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ABSTRACT 


This  research  investigated  multiple  detonation  diffraction  events  in  order  to  better 
understand  the  limits  and  benefits  of  diffraction  strategies  with  respect  to  pulse 
detonation  engine  design.  Hydrogen/air  detonations  were  generated  using  swept 
ramp  obstacles  in  a  1 .27  m  long  channel  with  a  cross  section  of  25.4  mm  by  88.9 
mm  and  were  diffracted  into  various  multiple-stepped  openings.  This  allowed  the 
detonation  wave  diffraction  transmission  limits  to  be  determined  for  hydrogen/air 
mixtures  and  to  better  understand  reinitiating  mechanisms  throughout  the 
diffraction  process.  Tests  were  conducted  for  area  ratios  ranging  from  2.00-2.60 
with  varying  equivalence  ratios  from  0.5-1 .5. 

Computational  methods  were  used  to  better  understand  the  diffraction 
phenomenon  using  a  series  of  sensitivity  studies  for  different  chemistry  sets, 
computational  cell  size  and  equivalence  ratio.  Experimental  tests  used  combined 
optical  shadowgraph  and  particle  image  velocimetry  imaging  systems  to  provide 
shock  wave  detail  and  velocity  information.  The  images  were  observed  through 
a  newly  designed  explosive  proof  optical  section  and  split  flow  detonation 
channel.  It  was  found  that  area  ratios  of  2.0  could  survive  single  and  double 
diffraction  events  over  a  range  an  equivalence  ratio  range  of  0.8  to  1.14  Area 
ratios  of  2.3  survived  the  primary  diffraction  event  for  equivalence  ratios  near 
stoichiometric  for  the  given  step  length.  Detonation  diffraction  for  area  ratios  of 
2.6  did  not  survive  the  primary  diffraction  event  for  any  equivalence  ratio  and 
were  unable  to  transmit  to  a  larger  combustor. 
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I.  INTRODUCTION 


It  is  estimated  that  nearly  80%  of  the  world's  energy  is  derived  from  the 
combustion  of  hydrocarbon  fuels  in  various  forms.  In  today's  energy  conscience 
and  profit-driven  society,  even  modest  improvements  in  current  systems  can 
yield  staggering  economic  gains.  The  quest  to  continuously  improve  energy 
efficiencies  of  such  combustion  driven  systems  has  spawned  research  in 
numerous  areas.  Pressure  gain  combustion  schemes  such  as  detonation  and 
constant  volume  combustion  hold  considerable  promise  in  achieving  said 
efficiencies  [1].  They  have  inherently  higher  theoretical  thermal  efficiencies 
primarily  due  to  their  near  constant  volume  detonation  based  combustion  cycle 
which  generate  much  less  entropy  when  compared  to  constant  pressure  events, 
creating  a  pressure  rise  during  heat  release.  Pulse  detonation  combustors 
(PDC)  exploit  such  efficiencies.  When  applied  in  a  pulse  detonation  engine  they 
are  estimated  to  achieve  greater  specific  impulses  over  many  conventional 
systems  for  a  greater  range  of  flight  regimes,  as  seen  in  Figure  1 . 


Figure  1 .  Comparison  of  high-speed  propulsion  technologies  (From  [2]) 


1 


Unfortunately,  the  ability  to  detonate  high  mass  flow  rates  in  a  reliable  and 
repeatable  manner  has  been  shown  to  be  a  tremendous  challenge  [3].  This  is 
especially  true  for  accepted  kerosene-based  aviation  fuels,  which  are  much  more 
difficult  to  detonate  than  hydrogen.  Additionally,  the  energy  required  to  directly 
initiate  a  detonation  in  a  combustor  large  enough  to  be  of  practical  use  in  a 
pulsed  detonation  engine  can  be  very  difficult.  Naturally,  much  research  has 
attempted  to  negate  complex,  high-energy  ignition  mechanisms  by  using  smaller 
channels  or  predetonators  followed  by  a  transition  into  a  larger  section. 
Predetonators  or  initiators  typically  involve  transitioning  a  deflagration  into  a 
detonation  through  turbulence  devices  to  accelerate  the  flame.  Such  devices 
can  have  losses,  so  once  the  detonation  is  formed,  it  is  normally  preferable  to 
transition  a  wave  to  a  larger  cross  sectional  area.  This  allows  for  additional  fuel 
processing  by  a  surviving  detonation  wave.  Numerous  geometries  and 
combinations  have  been  used  to  date,  but  it  remains  to  be  seen  whether  a 
specific  geometry  can  offer  any  advantages  over  others.  Hence,  the  requirement 
to  better  understand  the  detonation  diffraction  phenomena  all  types  of 
geometries. 

Many  configurations  investigated  in  the  open  literature  include  Brophy  et 
al.,  [4,  5],  who  studied  a  number  of  initiator-combustor  combinations,  shown  in 
Figure  2a  through  Figure  2c,  with  oxygen-based  initiators  with  some  success. 


Figure  2.  Various  initiator  main  combustor  conjurations  (From  [3]) 
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Oxygen-based  diffracting  initiators  were  eventually  abandoned  due  to  U.S.  Navy 
safety  concerns  and  ultimately  led  the  way  to  a  new  generation  of  swept  ramp 
initiators,  which  have  been  shown  to  reduce  pressure  losses  by  a  factor  of  5  over 
conventional  turbulence  devices  [6]. 

A  great  deal  of  historical  detonation  research  deals  with  fundamental 
diffraction  cases  into  unconfined  spaces  and  provides  tremendous  understanding 
of  the  detonation  phenomenon.  Unfortunately,  there  is  considerably  less 
literature  on  subcritical  reinitiating  events,  whereby  initiator  sections  are  made  as 
small  possible  to  support  a  multiheaded  detonation  wave  front.  Such  cases  are 
of  prime  concern  in  the  aero  propulsion  world,  since  weight  and  space  are  at  a 
premium.  In  order  to  diffract  detonations  from  heavier  and  more  stable  fuels  with 
larger  cell  structures  such  as  JP-10,  initiators  will  no  doubt  be  required  to 
operate  at  elevated  pressures  in  a  limited  amount  of  volume,  so  the  need  to 
understand  how  aggressive  a  detonation  can  diffract  is  very  important. 

This  research  aims  to  continue  investigations  in  detonation  diffraction  by 
exploring  multiple  diffractions  in  a  1  to  3.5  aspect  ratio  channel,  in  order  to  better 
understand  the  limits  and  benefits  of  diffraction  strategies  with  this  respect  to 
pulse  detonation  engine  design.  Hydrogen/air  detonations  were  generated  using 
swept  ramp  obstacles  in  a  1.27  m  long  channel  with  a  cross  section  of  25.4  mm 
by  88.9  mm  and  were  diffracted  into  various  multiple-stepped  openings.  Tests 
were  conducted  for  area  ratios  ranging  from  2.00-2.60,  with  varying  equivalence 
ratios  from  0.5-1 .5. 

Computational  methods  were  also  used  to  better  understand  the 
diffraction  phenomenon  using  a  series  of  sensitivity  studies  for  different  reduced 
chemistry  sets,  computational  cell  size  and  equivalence  ratio.  Experimental  tests 
used  combined  optical  shadowgraph  and  particle  image  velocimetry  imaging 
systems  to  provide  shock  wave  detail  and  velocity  information. 

Results  from  the  research  provide  baseline  design  criteria  for  some 
current  pulse  detonation  engine  concepts  being  investigated  in  Figure  3. 
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Figure  3.  Cross  sectional  view  of  a  Concentric  Tube  and  multiple  high  aspect 

ratio  channel  detonation  combustors 


The  use  of  concentric  tubes  or  multiple  high  aspect  ratio  channels  to  transition 
additional  fuel  to  the  combustor  while  isolating  the  initiator  section  may  allow  for 
greater  thrust  production.  In  each  diffraction  plane  the  detonation  front  is 
supported  by  additional  air  and  fuel  mixture.  Such  concepts  have  the  potential  to 
produce  higher  specific  impulse  values  while  also  smoothing  the  inherent 
intermittent  production  of  thrust.  For  the  concept  to  work,  the  detonation  front 
must  survive  every  diffraction  event  through  an  appropriate  area  change. 
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II.  BACKGROUND 


A.  COMBUSTION  WAVE  TYPES 

Combustion  events  can  be  classified  into  three  categories,  explosion, 
deflagration  and  detonations.  Generally,  in  explosions  the  rate  of  heat 
generation  is  extremely  fast  and  does  not  require  a  combustion  wave  through  the 
exploding  media.  Deflagrations  on  the  other  hand  propagate  at  much  lower 
velocities  with  regard  to  reactants  they  transmit  through.  They  are  essentially 
subsonic  expansions  waves  with  pressure  drops  across  the  reaction  zone  due  to 
products  that  are  accelerated  away  from  the  wave  in  the  opposite  direction 
proportional  to  diffusivity  and  rate  of  reaction.  Detonations  are  characterized  by 
supersonic  combustion  waves  that  propagate  extremely  fast  through  a  reactive 
mixture  increasing  the  thermodynamic  properties  sharply,  creating  a  shock- 
induced  chemical  reaction,  which  drives  the  detonation.  Some  basic  detonation 
models  include  Chapman-Jouguet  and  the  ZND  models.  Main  differences  are 
depicted  in  Table  1. 


Table  1.  Detonation  and  deflagration  parameters  (From  [1]) 


Detonation 

Deflagration 

ul/cl  (Mach  Number) 

5-10 

0.0001  -0.03 

u2/u1 

0.4  -  0.7  (deceleration) 

4-6  (acceleration) 

P2/P1 

13-55  (compression) 

*  0.98  (slight  expansion) 

T2/T1 

8-21  (heat  addition) 

4-16  (heat  addition) 

p2/p1 

1.7 -2.6 

0.06-0.25 

B.  RANKINE  HUGONIOT  GAS  DYANAMIC  RELATIONS 

Rankine  (1870)  and  Hugoniot  (1887-1889)  were  the  first  to  investigate 
solutions  to  the  non-reacting  conservation  equations  involving  shockwave,  often 
termed  the  Rankine-Hugoniot  jump  conditions  since  they  applied  a  control 
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volume  across  the  shock  looking  solely  at  the  upstream  and  downstream 
components.  In  doing  so,  they  were  able  to  derive  theses  components  in  terms 
of  wave  speed,  pressures  and  particle  velocity  post  shockwave.  Lee  [7]  provides 
a  nice  derivation  of  these  relations,  starting  with  the  one-dimensional  steady 
state  conservation  equations  across  a  fixed  detonation  wave,  specified  as  follows 
with  subscripts  0  and  1  representing  reactant  and  product  states  with  velocities 
determined  with  respect  to  the  wave,  as  depicted  in  Figure  4.  This  theory 
provides  a  good  first  approximation  to  the  changes  in  properties  during 
detonations  and  should  also  be  noted  that  these  equations  assume  a  steady  and 
adiabatic  front. 


Ui,Pi,Pi,h  1 


> 


Figure  4.  One-dimension  steady  state  model 


MaSS  jSg.'lig.  —  (1) 

Momentum  -t-  (2) 

Energy  q  +  ^  -  ti2  +  ^  (3) 

If  we  specify  an  equation  of  state  for  both  reactant  and  products  using  the  perfect 
gas  law  ( j?  =  pfir),  and  if  we  assume  &  =  with  the  relationship  —  Cv  =  R 

and  f  enthalpy  can  now  be  written  as  follows: 

A1 

Y  P 

r  -  ip  (4) 
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Since  we  have  five  unknowns  (wi/fvfv&i  «-■»££%),  the  latter  being  detonation 
wave  speed  and  only  four  equations,  one  more  equation  is  required  to  close  the 
set. 

In  order  to  close  this  equation,  we  need  to  define  the  Hugoniot  curve 
equations.  If  we  start  by  combining  Equations  1  and  2,  we  obtain: 


(5) 


(6) 


In  this  case,  v  =  l/p  is  known  as  specific  volume  and  rti  =  pit  mass  flux  per  unit 
area.  If  we  start  by  observing  the  pressure-specific  volume  plane  of  all 
theoretically  possible  solutions,  we  can  start  to  better  understand  relations 
between  Rayleigh  lines  and  Hugoniot  curves.  Looking  at  Equation  6,  we  can  see 
that  real  solutions  only  exist  in  the  detonation  and  deflagration  regions  in  white  in 
Figure  5. 
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Solutions  in  the  p-v  plane  (From  [7]) 


For  solutions  to  exist,  we  have  either  j?a  >  pa,Vj}  >vx  ,  the  case  where  we 
have  a  compression  solution  for  a  detonation  or  p1  an  expansion 

solution  and  deflagration,  as  in  Figure  5  Lee  defines  x  =  =  pQ/p±  and 

7  =  Pi/ft,  in  doing  so  one  can  now  rewrite  Equation  6  as: 

7  =  (1  -  VaMa)  ~  (7) 


this  is  the  Rayleigh  line  equation  with  a  slope  of  we  observe  that  the 

combustion  wave  velocity  is  directly  proportional  to  square  root  of  the  slope. 
Equation  7  describes  a  thermodynamic  path  the  gas  takes  to  across  the 
combustion  wave  from  state  (1 ,  1 )  to  state  (x,  y). 
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The  Hugoniot  curve  represents  an  agglomeration  of  all  possible  attainable 
downstream  states,  for  a  given  specific  set  of  initial  conditions,  as  depicted  in 
Figure  6.  By  eliminating  velocities  in  Equation  3,  we  obtain  the  following 
equation  for  the  Hugoniot  curve: 


1 

hx  -  +  $  -  j +  J*j) fa,  -  (g) 

We  also  know  that  solutions  to  the  conservation  equations  must  satisfy  both 
Rayleigh  and  Hugoniot  curves  simultaneously.  In  other  words,  the  evolution  of 
reactants  (x=y=1)  to  products  state  (x,  y)  follow  a  Rayleigh  line  until  intersection 
on  a  Hugoniot  curve,  seen  again  at  Figure  6.  For  finite  heat  addition  q,  the 
Rayleigh  line  will  intersect  four  points  on  the  Hugoniot  curve,  two  in  each  of  the 
detonation  and  deflagration  branches.  In  the  strong  detonation  region,  both 
pressure  and  density  increase  drastically  and  require  overdriven  initial  conditions. 
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Figure  6.  Rayleigh  line  and  Hugoniot  curve  (From  [7]) 


Two  other  solutions  exist  when  the  Rayleigh  line  is  tangent  to  the 
Hugoniot  curve.  In  this  case,  we  obtain  a  minimum  velocity  on  the  detonation 
branch  and  a  Mach  number  of  the  burned  gas  of  Mi=1.  The  tangency  solutions 
are  known  as  Chapman-Jouguet  solutions  seen  in  Figure  7.  The  upper  solution 
agrees  quite  well  to  the  experimental  data  and  is  easily  calculated  using  Newton- 
Raphson  iterations  and  has  been  incorporated  in  such  programs  as  the  NASA 
“Chemical  Equilibrium  with  Applications"  (CEA),  developed  by  McBride,  Zehe 
and  Gordon  [9]. 
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Figure  7.  Chapman-Jouguet  tangency  solutions  (From  [7]) 


C.  ZEL’DOVICH-VON  NEUMANN-DORING  (ZND)  ONE  DIMENSIONAL 

WAVE  STRUCTURE 

A  simple  and  widely  accepted  model  of  a  detonation  wave  is  the  ZND 
model  comprising  of  a  strong  shock  front  followed  by  a  much  thicker  induction 
and  reaction  zone  travelling  at  Chapman-Jouguet  velocity.  The  shockwave 
initially  adiabatically  compresses  the  mixture  to  a  peak  pressure  known  as  the 
Van  Neumann  state,  enabling  dissociation  of  the  molecules  within  a  thermally 
neutral  induction  zone.  At  a  certain  point,  enough  free  radicals  are  created  to 
start  the  chain  branching  reactions  eventually  converting  all  reactants  into 
products  within  a  reaction  zone.  The  rapid  chemical  reaction  rates  enable  a  very 
quick  rise  in  temperature  dropping  both  pressure  and  density  in  the  reaction  zone 
the  resulting  expansion  accelerate  the  gases  away  behind  the  front  coupling  it 
the  acoustic  front  [7]. 
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(a)  Steady,  one -dimensional  ZNT)  model. 


(b)  Laser  shadowgraph  of  self-propa¬ 
gating  detonation  wave 
iH2  +  0.5O2.  P]  =  20  kPa) 


Figure  8.  Model  and  Real  detonation  (From  [8]) 


D.  DETONATION  CELL  STRUCTURE 

Detonation  fronts  in  real  life  have  complex  three  dimensional  structures. 
The  wave  front  contains  multiple  and  constantly  interacting  curved  incident,  Mach 
stems  and  reflected  shock  waves  meet  at  so  called  triple  points.  These 
structures  are  highly  transient  due  to  spatial  and  temporal  instabilities.  If  we 
observe  the  highly  resolved  Schlieren  photograph  on  Figure  9,  we  can  see  what 
the  reaction  looks  like  before  two  transverse  waves  collide.  The  accompanying 
sketch  provides  some  additional  explanation  into  all  the  interconnected  events. 
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..Mach  shock 


Burned  gases 


High-density  unbumed  gasa: 
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_ _ Triple-  point 
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Mach  shock 


Figure  9.  Schlieren  photograph  of  a  detonation  in  a  thin  channel  (From  [7]) 

The  shock  front  has  multiple  colliding  waves  that  progress  transversely 
across  the  front,  creating  a  triple  point  path  that  forms  a  diamond  shape.  The 
chemical  reactions  are  not  always  perfect  during  a  detonation  event;  note  the 
turbulent  nature  of  the  reaction  zone  behind  the  detonation  with  numerous  shear 
layers  and  unburned  pockets  that  are  consumed.  In  highly  unstable  mixtures  like 
hydrogen-air  and  many  hydrocarbons  have  strong  transverse  waves  which  tend 
to  drive  the  detonation.  This  can  reduce  the  amount  of  chemical  energy 
contributing  to  the  main  front  leaving  unburned  pockets  to  be  consumed  later  by 
turbulent  diffusion.  These  mixtures  can  also  exhibit  irregular  cell  structures  and 
irregular  triple  point  trajectories.  Contrary  to  some  highly  regular  cell  structured 
mixtures  that  contain  heavy  amounts  of  argon  as  depicted  in  Figure  10. 
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Detonations  with  regular  cell  structure  are  less  reliant  on  transverse  wave 
collisions  to  continue  the  propagation  of  the  front.  Generally,  these  mixtures  are 
more  reactive. 


Triple  point  trajectories 


Mach  stem 

Chemical  induction  layer 
Incident  shock 


Shear  layer 


Transverse  shock 


Figure  10.  Example  of  a  regular  triple  point  cell  structure  of  a  detonation 

(From  [7]) 


Triple  point  trajectories  can  be  mapped  by  soot  foil  techniques  whereby 
shear  and  pressure  gradients  rearrange  soot  particles  leaving  behind  evidence  of 
triple  point  paths.  From  these  soot  foils  traces,  cell  width  A  can  be  determined 
and  provides  for  an  estimation  of  mixture  reactivity.  The  smaller  the  cell  width  is 
the  tighter  the  coupling  is  between  chemical  kinetics  and  gas  dynamic  effects.  It 
should  be  noted  however  that  there  is  a  large  variance  in  measurements  in  cell 
structure  via  soot  foils  as  observed  by  Shepherd  [10].  Soot  foil  measurements 
can  be  very  subjective  with  errors  up  to  50%. 
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E. 


STOICHIOMETRY 


Combustion  conditions  of  fuel-oxidizer  mixtures  are  generally  specified  in 
terms  of  equivalence  ratio  <t>.  This  characterizes  the  proportion  of  fuel  to  oxidizer 
in  the  mixture.  Defined  as  the  ratio  of  the  fuel-oxidizer  ratio  to  the  stoichiometric 
fuel-oxidizer  ratio  and  can  be  mathematically  represented  as  follows: 

^  _  m  fuel  j  m oxidizer  _  fuel  j  W oxidizer 

^  fue^  ■'  oxidizer  '}sr0jchj0metjil-  ^  fv&l  /  oxidizer  ^ stojc]]j0mefi-jc  (9) 

The  composition  and  thermodynamic  properties  of  fuel-oxidizer  mixtures  during 
combustion  and  detonation  are  dependent  on  the  equivalence  ratio. 
Consequently  the  detonation  cell  width  is  also  highly  dependent  on  mixture  fuel 
ratio.  A  ratio  greater  than  1  implies  excess  fuel  or  fuel  rich  mixtures  than  would 
be  required  for  a  stoichiometric  reaction  and  conversely  a  ratio  less  than  one 
implies  excess  oxidizer  or  fuel  lean  mixtures.  Looking  more  closely  at  the 
dependence  of  detonation  cell  width  by  inspecting  Figure  11,  we  can  see  a  very 
strong  relation  to  equivalence  ratio.  Cell  width  for  a  particular  mixture  is  smallest 
at  or  very  close  to  0=1 . 
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Figure  1 1 .  Variation  in  detonation  cell  size  with  fuel  air  composition  (From  [1]) 

F.  DETONATION  DIFFRACTION 

Diffraction  occurs  when  a  confined  detonation  is  allowed  to  expand  into  a 
larger  cross  sectional  area  at  which  point  expansion  waves  are  communicated  to 
the  flow  from  the  behind  the  front.  Diverging  transverse  waves  affected  by  this 
perturbation  produce  unsteadiness  and  curvature  in  the  front,  which  by  definition 
also  reduces  chemical  reaction  rates.  This  reduces  the  energy  is  release  rate 
thereby  providing  less  support  to  the  detonation  front.  Lee  illustrates  this  in 
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Figure  12,  we  see  a  rarefaction  fan  eventually  propagates  toward  the  axis  of  the 
channel  or  tube  and  the  flame  front  decouples  from  the  acoustic  front  in  a  total 
failure. 


Figure  12.  Rarefaction  fan  penetrating  propagating  a  toward  detonation  axis 

(From  [8]) 


Whether  a  detonation  survives  or  fails  an  expansion  depends  on  a  number  of 
factors:  mixture  composition,  thermodynamic  initial  conditions,  confining 
geometry  or  boundary  conditions  and  detonation  velocity  at  diffraction  plane. 
Three  different  cases  have  been  observed:  super  critical  where  a  detonation 
successfully  transitions  a  diffraction  plane,  subcritical  cases  occur  when  a 
detonation  wave  is  quenched  and  critical  cases  are  known  as  the  transition  point 
between  the  two  previous  cases,  seen  in  Figure  13. 
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(a)  Supercritical:  suc¬ 
cessful  transmission. 


(b)  Critical:  re- initiation 
as  shock  decouples 
from  reaction  zone. 


(c)  Subcritical:  complete 
failure  as  shock  decou¬ 
ples  from  reaction  zone. 


Figure  13.  Diffraction  regimes  (From  [8]) 


The  detonation  wave  diffraction  phenomena  is  not  new;  however,  the  last 
75-plus  years  of  research  in  this  area  has  still  has  not  yielded  a  complete 
understanding  of  all  the  fluid  dynamic  complexities  in  this  event.  With  all 
detonations,  high  pressures  and  temperatures,  short  and  varying  time  scales 
have  stymied  much  of  the  research  involved  in  this  area.  Lafitte  [11]  was  the  first 
to  attempt  the  transition  of  a  detonation  wave  from  a  7mm  tube  into  a  spherical 
chamber  using  a  mixture  of  CS2+C>2.  He  was  unsuccessful.  Zeldovich's  [12] 
Streak  camera  experiments  exhibited  detonations  that  decayed  completely  into 
flames  for  some  conditions  while  others  continued  to  detonate  across  the 
diffraction  plane,  he  established  that  for  various  tube  diameters  a  critical  diameter 
exists  for  successful  reinitiation  into  a  larger  space,  based  on  a  mixtures  initial 
conditions  and  confinement  properties.  With  improvements  in  experimental 
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techniques  in  the  1950s  through  1970s,  cellular  structure  of  the  detonation  wave 
was  better  characterized.  Shchelkin  and  Troshin  [13],  were  the  first  to  relate  cell 
size  to  the  Zeldovich-Newmann-Doring  (ZND)  model  reaction  length. 
Mitrofanov  and  Solouhkin  [14]  noted  the  first  empirical  correlation  relating 
detonation  cell  width  to  critical  tube  diameter;  they  proposed  that  the  critical 
diameter  is  13  times  the  detonation  cell  width  A.  They  demonstrated  that  for  a 
number  of  hydrocarbon  fuels  that  if  a  detonation  front  had  13  or  more  cell  widths 
across  a  tube  diameter,  the  detonation  would  survive  a  diffraction  event  into  any 
open  space  with  the  same  chemistry.  Follow-on  research  [15,  16,  17]  finds  this 
theory  breaks  down  for  more  regular  cell  structured  fuel  mixtures  and  higher 
aspect  ratio  rectangle  configurations  as  seen  in  Figure  14.  Extreme  aspect  ratio 
channels  great  than  5  asymptote  to  critical  channel  widths  Wc  of  4  times  the 
detonation  cell  width  A.  Such  channels  appear  to  more  easily  transition 
detonation  waves  as  other  transverse  modes  are  suppressed  and  expansion  is 
restricted  to  two  dimensions. 


Figure  14.  Critical  channel  height  as  a  function  of  aspect  ratio  (From  [7]) 
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Cell  structure  and  chemical  kinetics,  as  well  as  initial  and  boundary  condition, 
play  a  major  role  in  how  detonations  waves  expand,  break  down  and  reinitiate, 
as  well  as  the  initial  and  boundary  conditions  [17,  18]. 
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III.  COMPUTATIONAL  SETUP 


A.  INTRODUCTION 

Computational  fluid  dynamic  (CFD)  software  was  used  extensively  during 
the  research  in  order  to  compliment  experimental  results  and  to  better 
understand  all  the  complexities  of  a  diffracting  gaseous  detonation  wave. 
Numerical  results  also  provided  additional  details  of  the  flow  field  that  were 
difficult  or  impossible  to  do  so  experimentally.  Simulations  were  largely  based  on 
solving  the  two  dimensional  non-linear  compressible  inviscid  reactive 
conservation  equations  known  as  the  Reactive  Euler  equations  below  in 
conjunction  with  multi-step  reduced  chemistry  sets  with  nine  species  (H2,  O2,  H, 
O,  OH,  HO2,  H2O2,  H2O,  and  N2)  and  18/19  elementary  reactions  in  rectangular 
coordinates.. 

dp  d(pu)  a(pv)  _ 

3 1  +  dx  dy 

d(pu)  3  {p  +  pu2)  d(puv) 

3 1  dx  dy 
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In  short  conservation  of  mass,  momentum  in  both  x  and  y  directions,  energy 
along  with  all  species  equations  in  terms  of  Arrhenius  rates  and  an  equation  of 
state  were  used.  Metacomp’s  CFD++  software  was  used  to  solve  this  equation 
set.  Structured  meshes  were  generated  using  Pointwise  meshing  software 
package.  In  order  to  set-up  a  baseline  for  the  study,  a  number  of  sensitivity 
studies  were  conducted  using  computational  cell  size,  degree  of  overdrive, 
stoichiometry  and  chemistry.  Sensitivity  simulations  were  obtained  using  both 
implicit  and  explicit  schemes.  Due  to  the  transient  nature  of  detonations,  implicit 
simulations,  a  CFL  of  1  was  used  to  ensure  maximum  stability.  In  all  explicit 
time-marching  schemes,  CFL  number  was  ramped  from  0.5-1  during  the  first 
1500  iterations  allowing  the  simulation  to  resolve  initial  deflagration  to  detonation 
formation.  Consequently,  the  time  step  selection  in  explicit  time-marching 
simulations  were  critical  for  stability,  in  our  case  as  a  detonation  wave  crosses  a 
discrete  grid,  the  time  step  must  be  small  enough  to  resolve  it. 

B.  REDUCED  CHEMISTRY 

In  order  to  correctly  predict  where  and  how  detonation  waves  would 
diffract  and  whether  or  not  the  end  state  survived  or  failed,  heavy  reliance  on 
reduced  chemistry  was  required.  Simple  H2  and  02  reactions  systems  have 
been  known  to  encompass  hundreds  of  reactions.  However,  computational 
resources  are  normally  limited.  There  has  been  significant  research  in  using 
reduced  chemical  models  that  capture  the  most  important  characteristics  with 
respect  to  a  particular  application.  Models  in  the  past  have  either  assumed  one- 
step  equilibrium  based  chemistry  having  limited  success  in  predicting 
detonations.  This  has  tended  to  over  simplify  the  reaction  and  over  predict 
product  temperature.  Computational  power  today  has  allowed  for  multistep  finite 
rate  chemistry  allowing  better  prediction  of  heat  release  and  ignition  delay  times 
(induction  time)  during  a  detonation. 

This  research  looked  at  three  reduced  sets  that  have  less  than  20 
reactions,  none  of  which  were  developed  exclusively  for  detonation  based 
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combustion.  The  first  set  developed  by  Kurkov  &  Burrows  at  NASA  for  scramjets 
was  based  on  lower  combustion  pressure  regimes  [19].  The  second  set  was  a 
modified  Kurkov  set  that  included  the  some  energy  absorbing  NO  chemistry 
based  on  a  reduction  scheme  from  Kuo  [1].  The  last  set  used  was  developed  at 
Stanford  by  Peterson  and  Hanson  for  ram  accelerators  for  a  range  of  higher 
pressures  for  up  to  pressured  of  50  atm  [20],  as  seen  in  Table  2.  Hydrogen  air 
detonations  as  explained  earlier  are  highly  unstable  reactions  and  have 
tremendous  pressure  fluctuations  that  can  range  from  10  atm  to  27  atm  or  from 
0.5PCj  to  1.5PCj.  [7],  The  reaction  rates  were  expressed  in  Arrhenius  form  as 
follows: 


(18) 


The  reaction  rate  constant,  k  was  a  function  of  temperature,  pre-exponential 
parameter  A  described  the  frequency  of  collisions  for  the  species  that  result  in  a 
reaction  as  well  as  the  preferred  orientation  of  the  molecular  species  that  favor  a 
reaction.  The  exponent  n  expresses  the  temperature  dependence  of  the  pre¬ 
exponential  parameter.  Activation  energy,  Ea  defines  the  minimum  amount  of 
energy  required  for  a  reaction  to  occur  and  R  is  the  universal  gas  constant. 
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Table  2.  Peterson  &  Hanson  reduced  set  for  high  pressure  H2/02,  N2 
reaction  mechanism  for  ram  accelerator  (From  [20]) 


Rale  coefficient11 


Number 

Reaction3 

A 

n 

E 

Comments0 

1 

o+h2  =  h  +  oh 

5.00  x  104 

2.70 

6,290 

(R3) 

2 

h+o2  +  m=ho2+m 

2. SO  x  1018 

-0.90 

0 

(R33)d 

3 

H  4~  02  4-  02  =  HO  2  4-  Q2 

3.00  x  1020 

-1.70 

0 

(R34) 

4 

H  +  02  +  H20  =  HO?  +  H20 

9.38  x  10,s 

-0.80 

0 

(R35) 

5 

H  4-  02  +  N2  =  H0:  4-  N2 

2.60  x  10'* 

-1.20 

0 

(R36) 

6 

h+o2  =  o+oh 

8.30  x  1013 

0.00 

14413 

(R38) 

7 

H  +  H02  =  02  +  H2 

2.80  x  10'3 

0.00 

1,068 

(R45) 

8 

h  +  hg2  =  oh+oh 

1.34  x  1014 

0.00 

635 

(R46) 

9 

H  4-  H2  O?  =  H02  4-  H2 

1.21  x  107 

2.00 

5,200 

(R47) 

10 

oh  +  h2  =  h2o+h 

2.16  x  108 

1.50 

3,430 

(R84) 

11 

OH  +  OH  +  M  =  H202  +  M 

7.40  x  1 0'3 

-0.40 

0 

km f,  (R85)e  f 

2.30  x  10 18 

-0.90 

-1,700 

^0 

12 

OH  4-  H02  =0?  +H?0 

2.90  x  1013 

0.00 

-500 

(RS7) 

13 

OH-bH202=HO2  +  H2O 

1.75  x  10'2 

0.00 

320 

(R88JS 

5.80  x  1014 

0.00 

9,560 

kb,  (R89)s 

14 

HO?  4~  HO2  =  02  4“  H2O? 

1.30  x  10" 

0.00 

-1,630 

kc,  (R 1 1 5)h 

4.20  x  10'4 

0.00 

12.000 

*</,  (R1 16V1 

15 

q+g  +  a/=o2  +  m 

1.20  x  1G17 

-1.00 

0 

(Rl)1 

16 

0  +  H  +  M  =  0H  +  M 

5.00  x  1017 

-1.00 

0 

(R2)e 

17 

H  4-  OH  4-  M  =  H2  O  4-  M 

2.20  x  1022 

-2.00 

0 

(R43)i 

18 

h+h  +  m  =  h2  +  m 

1.00  x  10'8 

-1.00 

0 

(R39)k 

Note:  species  —  H,  O,  OH,  H2,  02,  H02.  H202,  H20,  and  N2  (or  Ar,  He,  etc. ). 

J  All  reactions  are  reversible. 

bfc(T)  =  AT"exp(— Ef  RT)\  units  are  in  cal,  mol,  cm3,  and  s. 

^  Corresponding  R  AMEC2*  reaction  numbers  are  in  parentheses. 

J M  does  not  include  02,  H20,  or  N2;  all  collision  efficiencies  =  1 .0. 

'Coll  is  ion  efficiencies  for  M\  N2  =  1.0,  Hi  =  2.0,  H20  =  6.0.  and  Ar  =  0.70, all  others  =  1.0. 

1  See  Table  2  for  True  parameters. 

pRate  coefficient  is  non-Arrhenius ;  Au  =  ka  +kb. 

h Rate  coefficient  is  non-Arrhenius;  frj 4  =  kc  +kj. 

’Collision  efficiencies  for  M :  N2  =  1 .0.  H2  =  2.4.  H20  =  1 5.4,  and  Ar  =  0.83,  all  others  =  1 .0. 
'Collision  efficiencies  for  Af:  N2  =  l.0,H2  =  0.73,  H20  =  3.65,  and  Ar  =  0.38,  all  others  =  1.0. 
k Co Ili sion  efficiencies  for  M:  N2  =  l.G.H?  =  1.7,  H20  =  7.0, and  Ar  =  0.63,  all  others  =  1.0. 


C.  COMPUTATIONAL  MESH 

The  computational  mesh  used  was  a  two  dimensional  channel  with  an 
open  step  in  the  positive  y  direction  with  a  1:6  height  expansion  ratio.  A 
resolution  of  50  pm  or  roughly  300  computation  domains  in  one  detonation  cell 
for  H2  was  used  equaling  to  approximately  4.5  million  cells.  The  mesh  was 
developed  using  Pointwise  64-bit  meshing  software  and  converted  to  run  in  a 
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CFD++  environment,  as  seen  in  Figure  15.  Some  researchers  have  reported 
that  the  resolution  requirement  to  fully  resolve  complex  chemistry  simulations  is 
an  element  size  of  only  0.1  pm  for  the  CJ  H2-air  detonations  at  1  atm  and  298  K 
[21].  Essentially  three  orders  of  magnitude  finer  of  the  induction  zone  thickness. 
Such  resolutions  would  be  well  beyond  currently  available  computational 
resources  and  could  not  be  attempted  in  this  study.  Additionally,  previous 
studies  at  NPS  [22]  have  shown  that  this  resolution  predicts  detonation  velocity 
within  5-10%  of  CJ  values.  It  was  hoped  that  this  resolution  would  be  sufficient 
to  capture  failure  modes  of  a  detonation  diffracting  around  an  open  step. 

The  model  was  scaled  down  slightly  from  the  original  experimental  setup, 
with  channel  height  scaled  down  to  1:2  and  the  channel  length  scaled  down  1:5, 
this  allowed  for  solutions  within  a  week. 
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D. 


BOUNDARY  AND  INITIAL  CONDITIONS 


Figure  16.  Boundary  and  Initial  Conditions  Computational  results 


Detonation  waves  in  the  simulations  were  created  by  introducing 
overdriven  high  pressure  products  in  small  region  at  the  beginning  into  a  main 
body  of  reactants.  Detonation  products  were  initiated  at  two  times  CJ  Pressure 
and  at  CJ  Temperature  conditions,  which  as  alluded  to  above  was  found  to  be 
the  minimum  overdrive  found  for  given  conditions  to  initiate  a  detonation  wave 
characteristic  of  the  fuel-oxidizer  mixture  as  part  of  the  initial  conditions.  Mass 
fraction  of  the  detonation  products  were  obtained  from  NASA  Chemical 
Equilibrium  Analysis  (CEA)  [9]  and  a  sample  is  included  in  Appendix  A.  Reactant 
initial  conditions  were  assumed  to  be  at  1  atm  and  298  K,  essentially  matching 
experimental  conditions  for  further  validation.  Typical  run  set  structure  was  as 
follows: 
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Equation  Set  Type 

•  Unsteady  Compressible  Euler 

•  Equation  of  State:  Ideal  Gas 


Riemann  Solver 

•  Minimum  Dissipation:  LHS 

•  No  pressure  switching 

Boundary  Conditions  (See  Figure  16) 

Time  Integration 

•  Implicit  CFL=1 

•  explicit  (for  some  runs)  based  a  CFL  ramped  from  0.5  -  1 
over  the  first  1500  iterations 

Spatial  Discretization 

•  2nd  Order  Accuracy  in  Space 

•  Dimension  of  polynomial:  2-D  XY 

•  Axis  stability  enhancement:  Yes 

•  Types  of  Total  Variation  Diminishing  (TVD)  Limiter:  Minmod 
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IV.  EXPERIMENTAL  SETUP 


A.  DETONATION  CHANNEL  MAIN  COMPONENTS 

The  combustor  design  was  driven  by  three  main  requirements,  the  desire 
to  use  hydrogen  and  air  as  reactants,  a  large  field  of  view  in  the  test  section  in 
order  to  facilitate  both  shadowgraph  and  particle  image  velocimetry  and  lastly, 
the  ability  to  adapt  the  existing  main  tube  design  in  test  cell  one.  The  main 
components  of  the  detonation  diffraction  test  setup  consists  of  an  inlet  manifold, 
a  48-inch  long  rectangular  channel  in  which  an  inner  split  flow  channel  is  inserted 
with  a  one  3.5-inch  tall  by  1-inch  wide  inner  channel  dimensions  and  an  optical 
section.  This  configuration  allows  for  detonations  to  pass  through  a  lower  section 
and  diffract  through  an  upper  section  carrying  unreacted  fuel-air  mixture  into  a 
viewable  optical  section  seen  at  Figure  17.  Appendix  B  includes  all  engineering 
drawings  for  the  split  flow  insert  and  optical  section  components  as  well  as  the 
new  manifold 
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B.  COMBUSTOR  OPERATION 

This  pulse  detonation  combustor  setup  operates  on  a  relatively  simple 
single  cycle  mode  in  order  to  facilitate  optical  experimental  observations  with 
shadowgraph  and  Particle  Image  Velocimetry  (PIV).  It  was  operated  on  an 
optics  table  consisting  of  two  electro  pneumatic  ball  valves  that  controlled  fuel 
and  seeding  air  along  with  a  main  air  duct  regulated  by  a  computer  controlled 
pressure  regulator.  Also,  various  check  valves,  two  hand  pressure  regulators  to 
set  precise  pressures  for  fuel  and  seeding  along  with  two  mass  metering  chokes. 
Two  lasers  were  employed  to  illuminate  the  detonation  front  for  image  capture  for 
their  respective  camera  system  using  a  continuous  Argon-Ion  for  shadowgraph 
and  a  double  pulsed  Nd-Yag  for  PIV,  seen  in  Figure  19.  Two  kistler  high-speed 
pressure  transducers  acquired  high  speed  data  and  were  used  to  determine 
detonation  velocity  as  well  as  trigger  the  PIV  system  via  a  BNC  555  pulse  delay 
generator.  Labview  commercial  control  and  automation  software  was  used  to 
operate  the  test  engine.  The  standard  sequence  consisted  of  filling  the  upper  on 
lower  channel  with  a  seeded  titanium  dioxide  hydrogen/air  mixture  and  igniting  it 
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and  subsequently  measuring  pressure,  detonation  velocity,  acquiring  high  speed 
shadowgraph  images  and  particle  image  velocimetry  images.  See  Appendix  C 
for  standard  operating  procedures. 


Figure  19.  Overall  experimental  setup 
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Figure  20.  Overall  experimental  Setup 


C.  OPTICAL  SECTION 

The  tube  was  attached  to  an  8.5-inch  tall  by  12.5-inch  long  and  4.5-inch 
wide  optical  section  with  opposing  5-inch  by  9-inch  viewing  areas  using  45- 
degree  chamfered  flush  fit  BK7  optical  glass.  Various  step  configurations  can  be 
installed  into  the  section.  The  design  was  based  on  a  previous  section  by 
Fludovich  [4],  however  the  new  section  was  designed  to  be  more  modular  with  a 
larger  optical  viewing  area  bolting  directly  to  the  detonation  channel  with  access 
ports  above  and  below,  seen  in  Figure  21.  The  entire  section  was  machined  into 
mainly  one  inch  thick  304  stainless  steel  sections  and  was  assembled  with 
socket  cap  head  screws.  The  strong  design  is  required  in  order  to  withstand  any 
bending  that  could  compromise  the  1.75  inch  BK7  optical  glass  it  was  designed 
to  encase. 
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Figure  21 .  Stainless  Steel  304  optical  section 


Finite  element  analysis  was  performed  on  the  BK7  glass  components  of 
the  older  and  newer  designs  to  ascertain  maximum  stress  and  deflection  under 
shock  loading  conditions  in  order  to  ensure  that  the  new  design  would  perform 
equally  or  better.  A  glass  thickness  of  4.45  cm  (1.75  in)  was  selected  to 
accommodate  the  detonation  impulse  pressures  of  up  to  36  atm  for  up  to  1  ps 
shock  loading.  Results  of  these  analyses  are  presented  in  Figures  22  and  23. 
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Figure  22.  Deflection  results  for  BK7  glass 
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Figure  23.  Principal  stress  results  for  the  BK7  glass 
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D.  LOWER  CHANNEL  SPLIT  FLOW  INSERT 

In  order  for  a  detonation  to  diffract  into  an  open  step  using  the  existing 
main  square  channel,  a  split  flow  channel  assembly  was  designed  to  easily  slide 
into  the  existing  hardware,  essentially  separating  the  flow  into  two  channels,  an 
upper  and  lower  allow  for  a  detonation  to  initiate.  The  detonations  were  initiated 
in  the  lower  section  expanding  into  an  adjacent  top  channel  at  the  diffraction 
plane  directly  observed  through  the  optical  section.  The  flame  was  transitioned 
into  a  detonation  using  10  interspaced  30-degree  swept  ramps  with  a  total 
blockage  ratio  of  40%.  This  section  was  machined  from  6061  Al  stock  and  is 
socket  cap  screw  assembled,  more  detailed  drawings  can  be  viewed  in 
Appendix  B. 


Figure  24.  Split  flow  insert. 
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E. 


IGNITION  SYSTEM 


The  detonation  is  initiated  by  a  Unison  Vision-50  variable  capacitative 
discharge  using  an  aviation  grade  spark  plug  mounted  at  the  head  end  of  the 
tube.  Prior  to  testing,  the  igniter  was  configured  manually  to  provide  a  2  Joule 
spark  and  was  remotely  triggered  by  the  Labview  software  (Figure  25). 


F.  FUEL  AND  AIR  DELIVERY  SYSTEM 

High  pressure  air  was  provided  to  the  test  cell  from  the  rocket  lab  facility 

air  system  with  hydrogen  provided  from  a  single  bottle.  Main  control  room 

computers  controlled  supply  pressures  for  both  air  and  hydrogen  with  the 

Tescom  ER3000  Version  2.0  software,  allowing  reactant  pressures  to  be  set 

remotely  as  required.  Fuel,  air  and  seed  delivery  revolved  around  the  main  six 

port  inlet  manifold  which  was  designed  to  adapt  to  the  existing  pair  of  1  inch  main 

air  lines.  Further  upstream  air  is  fed  through  a  1.5-inch  pipe  using  a  0.235-inch 

flow  metered  choke.  The  manifold  design  delivers  a  high  enough  mass  flow, 
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ensuring  swept  ramp  effectiveness  and  thus  reliable  detonations.  The  design 
allows  for  low  pressures  seed  particulate  penetration  into  the  main  flow  and 
ensured  the  seed  reservoir's  2  bar  pressure  limit  was  not  exceeded.  Main  air 
lines  were  broken  down  into  two  successive  pairs  of  half  inch  lines  where  fuel 
and  seed  was  fed  to  the  main  flow  in  two  stages  into  the  six-holed  manifold  that 
accepted  up  to  five  lines  with  one  hole  tapped  for  a  standard  turbine  plug.,  as 
depicted  in  Figure  26.  Hydrogen  was  supplied  using  half-inch  tubing  that  in  turn 
fed  two  main  air  arms  just  prior  to  the  four  line  split. 


Figure  26.  Inlet  manifold  and  feed  lines 


Unfortunately  hydrogen  and  seeding  could  not  be  regulated  using  the 
ER3000  systems  since  operating  pressures  were  lower  then  what  this  system 
was  capable  of  delivering  to  within  reasonable  errors.  Hand  ball  valve  regulators 
were  used  instead,  since  typical  operating  pressures  were  under  100  psig  and  10 


37 


psig,  respectively.  All  pneumatic  ball  valves  were  controlled  by  LabVIEW 
software  via  Crydom  control  solenoid  switches  located  in  an  electronics  cabinet 
located  in  the  test  cell  as  seen  in  Figures  27  and  28.  Choke  diameters  used 
during  this  testing  were  0.236  inch  for  the  air  supply  and  0.157  inch  for  the 
hydrogen  supply.  Two  check  valves  located  in  the  fuel  and  seed  lines  were 
installed  to  prevent  the  possibility  of  backflow  due  to  the  detonation  event. 


LabVIEW 


ER3000 


PXI-1000B 


~7 

. . ■■■■■ 

p 

m 

Nl 

H2  Pressure 

HP  Air  Transducer 

Combustible  Mixture  ®  Regulator 
Shop  Air  [jT]  Ball  Valve 

Control  Line  [^J  Choked  Orifice 

Data  Line  li=M  One-Wav  Valve 


Figure  27.  Fuel  and  air  delivery  system  (From  [23]) 


G.  INSTRUMENTATION 

Instrumentation  employed  all  through  the  testing  controlled  initiation  as 
well  as  data  collection.  Each  of  the  pressure  signals  were  routed  through  one  of 
four  National  Instruments  (Nl)  14-bit  PXI-6115  cards  mounted  in  the  Nl  PXI- 
1000B  chassis,  shown  in  the  upper  half  of  Figure  28.  This  chassis  interfaced 
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with  the  computers  in  the  control  room  through  the  Nl  PXI-MXI-4  PXI  Bridge  and 
was  capable  of  collecting  either  real-time  or  high-speed  buffered  data.  Optical 
data  imaged  from  the  high-speed  camera  and  PIV  systems  were  routed  to  a 
dedicated  desktop  computer  in  the  control  room. 


Figure  28.  National  Instruments  PXI-1000B  Chassis  (upper  section)  and 
Crydom  Control  Solenoid  Switches  (lower  section)  (From  [23]). 

H.  DYNAMIC  PRESSURE  TRANSDUCERS 

Inside  the  channel,  two  high-frequency  Kistler  dynamic  pressure 

transducers  were  mounted  as  a  pair  in  order  to  calculate  the  wave  velocity  based 

on  the  pressure  transients.  The  transducers  were  connected  to  Kistler  Type 

5010  Dual  Mode  amplifiers  which  route  the  data  to  two  Nl  PXI-6115  data  cards 

in  the  PXI-1000B  for  high-speed  data  collection.  Once  testing  was  initiated,  the 

cards  were  configured  to  begin  collecting  data  with  the  triggering  of  the 
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capacitive  discharge  igniter  and  collected  data  at  a  rate  of  500  kHz.  Data  was 
stored  via  card  buffer  and  later  saved  to  computers  located  in  the  control  room. 
Pressure  data  supplied  from  Kistler  transducers  provided  high-resolution 
confirmation  of  shock  passage. 

I.  CONTROL 

Every  test  sequence  was  initiated  and  controlled  using  Labview  software. 
All  pressure  transducers  were  linked  into  to  the  code  to  allow  for  active 
monitoring  during  every  test  sequence.  A  sample  control  panel  can  be  seen  in 
Figure  29. 


Test  Cell  1  Live 

Filename  and  path 
%  C:\JuilletData\lded0.txt 
Highspeed  file  path 
l  C:\JuilletData\HighSpeedldeclO.txt 


Start  Sequence 


Ignition 


& 


Kistlerl 

Kistler2 


‘+;B*i  ] 


Sample  Rate  (hz)  H2Fuel 

Pressure  (psig) 

1000 

0 

Samples  Main  Air  PreChoke 

20  Pressure  (psig) 


Min  value  (volts) 
o 


o 

Node4  Seeder  Air 
Pressure  (psig) 


Max  value  (volts)  0 


Sample  High  Speed  Rate  (hz) 
500000 

Samples  Highspeed 
200000 

Min  value  HS  (volts) 

-10 

Max  value  HS  (volts) 

10 


Figure  29.  Sample  Labview  control  panel 
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J.  SHADOWGRAPH 

Opposing  optical  windows  allow  for  the  transmission  of  a  shadowgraph 
light  source  to  pass  through  the  section  during  a  detonation.  Shadowgraph  setup 
utilizes  an  Ion  Laser  Technology  Model  55001  Argon  Ion  laser  emitting  750  mW 
beam  at  488  nm  passing  into  a  6.3  pm  spatial  filter  placed  in  front  of  the  laser 
aperture.  This  produced  a  spatially  filtered  light  source,  which  became 
collimated  after  being  reflected  from  a  concave  mirror.  The  resulting  beam  was 
directed  through  the  optical  test  section  using  a  flat  mirror.  On  the  opposite  side 
of  the  optical  test  section,  a  flat  mirror  directs  laser  light  toward  a  focusing  mirror 
placed  which  in  turn  is  directs  the  light  directly  into  the  ultra  17  intensified  high¬ 
speed  camera. 

The  main  rationale  for  choosing  shadowgraph  as  opposed  to  Schleiren 
technique  was  due  to  the  high  coherence  length  of  the  argon  ion  laser  system. 
Such  high  coherence  length,  do  not  provide  for  good  Schleiren  images.  The 
argon  ion  system  was  used  as  monochromatic  light  source. 

Invisible  Vision's  Ultra  17  high-speed  camera  suite  was  used  to  capture  all 
shadowgraph  images.  The  system  was  triggered  from  the  PI  kistler  high  speed 
pressure  probe  with  delays  ranging  from  of  225-275  ps  depending  on  what 
portion  of  the  optical  window  was  being  captured.  At  the  150,000  fps  mode,  a 
resolution  of  512  X  512  pixels  was  possible  with  this  system.  Exposure  settings 
varied  between  50-500  ns,  depending  on  the  amount  of  light  in  the  test  cell  while 
allowing  gain  to  be  lower  than  30%.  Inter  frame  times  were  6.17  ps. 
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Figure  30.  Shadowgraph  setup 


K.  PARTICLE  IMAGE  VELOCIMETRY  (PIV) 

The  PIV  arrangement  used  a  Litron  Nd:YAG  Dual  Cavity  pulsed  laser  2  x 
200  mJ  at  532  nm  with  up  to  15  Hz  pulse  rate  with  variable  light  optics.  The  laser 
sheet  was  passed  underneath  perpendicular  to  the  flow  through  a  narrow  optical 
slot  Image  as  seen  in  Appendix  B.  Image  capturing  was  done  using  a  Lavsion 
Imager  Pro  X  2M  1648  x  1214  pixel  resolution  camera  fitted  with  a  50  mm  f  1 .4 
Nikon  lens  and  532  nm  filter.  Laser  and  camera  timing  was  accomplished  using 
a  PC  based  programmable  timing  unit  (PTU-9).  All  image  post-processing  was 
completed  with  DaVis  software  installed  on  a  standard  quad  core  PC. 

In  order  to  properly  seed,  a  reacting  and  detonating  proper  seed  material 
would  be  required  to  withstand  high  temperature  ad  pressures  while  also 
exhibiting  minimum  interference  to  the  flow.  Micron  particle  Titanium  Dioxide 
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was  selected  for  its  good  stability  at  high  temperatures  and  its  relative  inertness. 
As  eluded  to  earlier  manifold  flow  seeding  was  accomplished  using  the  Lavision’s 
Particle  blaster  100  and  titanium  dioxide  seeding  particles.  Since  testing 
involved  simultaneous  PIV  and  shadowgraph,  the  PIV  camera  was  angled  out  of 
the  shadowgraph  optical  path,  requiring  calibration  within  the  DaVis  software 


Figure  31.  PIV  setup 


L.  OPTICAL  SECTION  TEST  CONFIGURATIONS 

There  were  five  major  test  configurations:  Configuration  1  consisted  of 
open  diffraction  with  no  step  sections  and  a  series  of  cold  flow  tests. 
Configuration  2,  3  and  4  were  single  step  diffraction  setups  of  various  height 
ratios;  5  and  6  were  double  step  diffraction  setups  with  two  different  height  ratios. 
Figure  32  depicts  common  setup  configurations. 
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Figure  32.  Single  and  double  step  configurations 


44 


V.  COMPUTATIONAL  RESULTS 


A.  OVERVIEW  OF  RESULTS 

It  was  originally  envisioned  that  the  CFD  portion  of  this  thesis  could 
potentially  shed  light  and  perhaps  help  us  better  understand  the  complexities  of 
gaseous  diffraction  detonation  wave  phenomena  while  also  assisting  in  the 
design  of  some  of  the  components.  Simulations  in  general  were  effective  in 
forming  a  detonation  given  the  various  parameters  but  once  the  detonation  front 
was  formed,  failure  of  the  front  during  diffraction  was  not  attained.  Unfortunately, 
not  one  simulation  was  observed  to  completely  fail  across  the  diffraction  plane 
and  thus  precluded  us  from  any  useful  results  in  computational  diffraction 
detonation. 

Two  main  reasons  were  suspected  for  simulations  that  could  not  resolve  a 
failing  detonation  front,  namely  poor  mesh  resolution  and  the  reduced  chemistry 
sets  utilized.  Poor  Resolution  (causing  artificial  reactivity)  resulted  in  an 
unresolved  the  triple  shock  system  at  the  front.  Incorrect  chemistry  could  not 
predict  the  proper  heat  release  and  ignition  delay  time  (induction  time),  as  well  as 
not  enough  diffusive  elements  in  the  simulation.  This  resulted  in  solutions  that 
were  overly  reactive  and  not  capable  of  failure  during  a  diffraction  event. 

Simulations  for  the  smallest  mesh  size  of  0.05  mm  used  for  the  study  ran 
for  roughly  240  hours  on  128  processors  on  the  new  NPS  hamming  Linux 
cluster,  while  the  largest  mesh  sizes  took  only  24  hours  to  run  with  the  same 
computational  settings.  CFL  number  for  all  simulations  was  in  between  0.5  and 
1  and  with  varying  time  steps  depending  on  the  mesh  size.  Table  3  is  a 
summary  of  all  the  computed  cases. 
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Table  3. 


Summary  of  computed  results  with  ATM  initial  conditions 


Run  # 

Mes 

h 

Size 

(mm) 

$ 

Reduced 
Chemistry  Set 

Numerical  Scheme  &  other 
settings 

Result 

1  — 23jul 

0.05 

1 

Kurkov  Burrows 

Implicit,  CFL  1 

No  failure 

2  -26jul 

0.1 

1 

Kurkov  Burrows 

Implicit,  CFL  1 

No  failure 

3- 

16aug 

0.1 

1 

Kurkov  Burrows 

Explicit,  CFL  ramped  0.5-1 

No  failure 

4- 

20aug 

0.05 

1 

Kurkov  Burrows 

Explicit,  CFL  1,  1C  1.5Pcj 

No 

detonation 

5- 

23aug 

0.05 

1 

Kurkov  Burrows 

Explicit,  CFL  1,  1C  2Pcj 

No  failure 

5- 

23aug 

0.05 

1 

Kurkov  Burrows 

Explicit,  CFL  1,  1C  2Pcj 

No  failure 

6- 

27aug 

0.05 

1 

Kurkov  Burrows 

Implicit,  CFL  1,  LHS 

dissipation 

No  failure 

7- 

30aug 

0.05 

1 

Kurkov  Burrows 

Implicit,  CFL  1,  LHS 
dissipation,  1C  1.5PCj, 

viscosity  on  walls 

No 

detonation 

8- 

31aug 

0.05 

1 

Kurkov,  Burrows 

Implicit,  CFL  1,  LHS 

dissipation,  1C  2PCj,  viscosity 
on  walls 

No  failure 

9-2sep 

0.05 

1 

Modified  Kurkov 
Burrows 

Implicit,  CFL  1,  LHS 

dissipation,  1C  2PCj,  added 
NO  chemistry 

No  failure 

10- 

7sep 

0.05 

1 

Modified  Kurkov 
Burrows 

Implicit,  CFL  1,  LHS 

dissipation,  1C  2.5PCj, 

viscosity  on  walls,  added  NO 
chemistry 

No  failure 

11- 

lOsep 

0.05 

0.8 

Modified  Kurkov 
Burrows 

Implicit,  CFL  1,  LHS 

dissipation,  1C  2PCj,  added 
NO  chemistry 

No  failure 

12- 

loct 

0.05 

0.8 

Modified  Kurkov 
Burrows 

Implicit,  CFL  1,  LHS 

dissipation,  1C  2PCj,  added 
NO  chemistry,  symmetry 
plane 

No  failure 

13- 

2oct 

0.05 

0.6 

Modified  Kurkov 
Burrows 

Implicit,  CFL  1,  LHS 

dissipation,  1C  2PCj,  added 
NO  chemistry,  symmetry 
plane 

No  failure 

14- 

12oct 

0.05 

0.6 

Modified  Kurkov 
Burrows 

Implicit,  CFL  1,  LHS 

dissipation,  1C  2PCj,  added 

No  failure 
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NO  chemistry 

15- 

0.05 

1 

Modified  Kurkov 

Explicit,  CFL  ramped  from  0.5 

No  failure 

26oct 

Burrows 

-0.95,  ran  faster,  2PCj,  added 
NO  chemistry 

16- 

0.05 

1 

Hanson, 

Explicit,  CFL  ramped  from  0.5 

No  failure 

5nov 

Peterson 

-0.95,  ran  faster,  2Pcj,  added 
NO  chemistry 

17- 

16nov 

0.05 

1 

Hanson, 

Peterson 

Explicit,  CFL  ramped  from  0.5 
-0.95,  ran  faster,  2PCj,  added 
NO  chemistry,  symmetry 
plane 

No  failure 

B.  DETONATION  FRONT  STRUCTURE 

The  following  figures  show  the  computed  contours  of  pressure,  OH  mass 
fraction,  pressure  and  temperature  plots  for  0.1  mm  and  0.05  mm  mesh  sizes. 
The  combustion  front  at  these  resolutions  did  not  capture  enough  front  details. 
Best  resolved  fronts  came  from  the  OH  contours  at  the  finest  resolution. 

Upon  first  inspection  at  the  detonation  front  on  the  OH  mass  fraction 
contours  right  before  the  diffraction  plane,  the  front  appears  to  lack  real 
sharpness  of  a  triple  point  that  we  would  normally  from  a  two  dimensional 
simulation.  Very  faint  keystone  structures  were  only  observed  with  a  0.05  mm 
mesh.  There  was  no  defined  zone  where  the  three  shocks  meet.  Looking  at 
Pintgen’s  extensive  PHD  work  using  planar  laser  induced  florescence  (PLIF)  of 
the  OH  radicals  he  revealed  their  presence  near  the  front  in  the  distinct  keystone 
patterns  is  clear  in  Figure  33. 
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Figure  33.  Observations  in  the  sub-critical  regime,  combined  Shlieren,  OH* 
and  PLIF  imaging  of  the  highly  irregular  H2-N2O  mixture  (From  [18]) 


The  H2  N2O  mixture  in  Pintgen  work  has  a  much  more  irregular  cell 
structure  and  is  much  more  reactive  than  our  H2  air  case  with  cell  sizes 
averaging  1-4  mm.  We  use  the  case  only  to  compare  irregular  and  subcritical 
type  trends.  Normally  during  high  temperature  and  pressure  events,  dissociation 
occurs  immediately  after  shock  passage  resulting  in  high  mass  fractions  in 
vicinity  of  the  front,  this  is  just  barely  apparent  with  these  contours.  A  run  with  a 
5pm  mesh  size  or  smaller  would  have  been  attempted  except  that  but  such  runs 
would  have  taken  months  for  results  with  NPS'  current  computational  resources 
and  would  not  have  been  feasible  given  current  time  constraints.  It  is  understood 
that  for  structured  meshes  a  5  pm  resolution  has  been  reported  as  the  starting 
point  for  adequate  detonation  front  structure  [21].  Looking  at  Figure  34  OH  mass 
fraction  contour  plots,  there  are  drastic  improvements  in  front  details  moving  from 
100  pm  to  a  50  pm  structured  mesh — regrettably,  not  enough  to  resolve  any 
jagged  edges  or  keystones  that  we  can  clearly  see  in  Pintgen's  results. 
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Figure  34.  Mass  fraction  OH  contour  plots  of  detonation  front  using  a  100  pm 

and  50  pm  structured  mesh 


With  the  pressure  contours  for  all  three  mesh  sizes,  we  can  measure 
detonation  cell  size  A  from  remnant  transverse  pressure  waves  right  behind  the 
front.  It  is  interesting  to  note  that  average  cell  size  in  each  mesh  size  seem  to  be 
consistently  smaller  then  experimental  values  of  1 .5  cm  for  Hydrogen.  For  our 
baseline  case  with  a  mesh  size  of  50  pm,  we  found  an  average  cell  size  of 
0.22  cm,  indicating  a  much  more  reactive  mixture.  Many  researchers  have 
reported  that  in  order  to  computationally  acquire  accurate  cell  size,  proper  grid 
resolution  is  critical.  Poor  resolution  has  also  been  known  to  cause  smaller  than 
normal  detonation  cell  sizes  and  in  other  words  a  much  more  reactive  mixture. 
Given  the  computational  constraints  we  could  not  confirm  this  since  we  did  not 
run  the  5pm  mesh  case.  Transverse  waves  also  revealed  some  additional 
irregularities;  they  seemed  to  persist  much  longer  after  the  front  had  passed  than 
normal.  They  occurred  much  more  regularly  then  one  would  expect  for  such  an 
irregular  mixture  like  H2-Air.  Pressures  in  and  around  the  front  were  on  average 
20%  higher  for  both  Von  Neumann  and  CJ  pressures.  This  was  expected  since 
nitrogen  chemistry  was  not  included. 
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Figure  35.  Pressure  contour  plots  of  detonation  front  using  a  100  pm  and 

50  pm  structured  mesh 


Synonymous  with  detonations  was  also  the  turbulent  nature  of  the 
reaction  zone.  Viscosity  and  transport  effects  played  an  extremely  vital  function 
within  the  small  scale  phenomena  of  the  detonation  structure.  None  of  this  detail 
is  apparent  at  this  mesh  size  in  any  of  the  contours.  No  shear  layers  were 
observable.  Metacomp,  maker  of  CFD++  stress  that  their  Euler  solver  settings 
employ  a  behind  the  scenes  dissipation  parameter  built  in  the  code.  We  could 
not  verify  whether  this  had  any  effect. 


Figure  36.  Temperature  contour  plots  of  detonation  front  using  a  100  pm  and 

50  pm  structured  mesh 
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Inspecting  H2  mass  fraction  contours  also  show  scarce  signs  of  unburned 
pockets  of  hydrogen  right  behind  the  front.  Based  on  the  literature,  this  has  been 
observed  experimentally  and  has  been  confirmed  computationally.  Lee  explains 
that  local  fluctuations  of  the  main  shock  front  during  the  low  velocity  phase  of  the 
cycle  can  increase  induction  times  resulting  in  later  heat  release  after  the  front. 
Many  irregular  fuels  can  have  significant  levels  of  unburnt  reactant. 

C.  CHEMICAL  KINETICS 

Initial  runs  used  the  Kurkov  and  Burrows  chemistry  set  for  hydrogen  as  an 
initial  investigation  to  see  whether  a  chemistry  set  developed  for  supersonic  low 
pressure  combustion  could  be  utilized  in  a  detonation  environment.  However,  we 
soon  discovered  this  set  was  not  the  most  ideal  set,  as  diffractions  set  forth  into 
unconfined  areas  would  not  completely  fail  with  our  smallest  structured  mesh. 
Granted  it  is  also  suspected  that  this  may  be  a  product  of  an  over  reactive 
mixture  at  the  current  grid  resolution.  Further  investigation  was  deemed  required 
at  this  point. 

The  Hanson  and  Peterson  set  was  the  last  iteration,  this  reduced  set  was 
exclusively  created  for  high  pressure  detonation  for  ram  accelerator  type 
combustion  events.  Unfortunately,  results  were  not  as  promising  as  hoped,  since 
an  average  was  taken  for  4  reaction  equations,  and  since  they  were  non- 
Arrhenius  and  could  not  be  readily  inputted  into  CFD++.  These  simulations  also 
were  unsuccessful  in  producing  a  failure  in  the  diffraction  front. 

Detonation  velocity  between  the  various  chemistry  sets  was  compared 
and  determined  for  a  number  of  runs  based  on  a  finite  difference  method  by 
measuring  the  front  change  in  position  between  two  time  steps,  see  Table  4 
some  the  typical  values  calculated.  Values  were  on  the  order  of  5%-10% 
percent  higher  than  the  CEA  calculated  values.  This  was  expected  as  the 
nitrogen  chemistry  was  omitted  for  all  Kurkov  &  Burrows  chemistry  sets.  Our 
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modified  Kurkov  set  included  a  one  step  NO  equation  which  was  hoped  to  slow 
down  the  chemical  kinetics.  Velocities  were  slower  but  the  code  still  was  not 
able  to  resolve  a  failing  diffraction. 


Table  4.  Calculated  detonation  velocity  with  CEA  calculated  values 


Run  # 

Mesh 

$ 

Reduced 

Numerical  Scheme  &  other 

Avg 

Size 

(mm) 

Chemistry 

Set 

settings 

Detonation 

Velocity 

(m/s) 

6- 

27aug 

0.05 

1 

Kurkov 

Burrows 

Implicit,  CFL  1,  LHS  dissipation 

2010 

7-30 

0.05 

1 

Kurkov 

Implicit,  CFL  1,  LHS 

2000 

Aug 

Burrows 

dissipation,  1C  1.5PCj,  viscosity 
on  walls 

11- 

0.05 

1 

Modified 

Implicit,  CFL  1,  LHS 

1910 

lOsep 

Kurkov 

Burrows 

dissipation,  1C  2Pcj,  added  NO 
chemistry 

CEA 

1 

NASA  Chemical  Equilibrium  VCj 

1964 

CEA 

0.8 

NASA  Chemical  Equilibrium  VCj 

1816 

D.  DETONATION  DIFFRACTION  CONTOUR  PLOTS 

All  simulations  exhibited  critical  diffraction  characteristics  even  though  all 
conditions  were  sub  critical.  All  detonations  undergo  some  partial  failure  which 
starts  almost  immediately  at  the  open  step  but  inevitably  a  re-ignition  event 
occurred  at  some  point,  reenergizing  the  entire  detonation.  Pressure, 
Temperature  and  OH  mass  fraction  contour  plots  at  one  mid  time  step  attest  to 
the  some  of  the  true  complexity  of  the  flow  system.  Near  the  main  axis  of  the 
channel  the  flame  front  never  really  appeared  to  decouple  from  the  precursor 
shock.  However,  the  flame  front  at  the  corner  completely  decouples.  Front 
speeds  were  observed  as  low  as  0.8VCj  near  the  corner  to  as  high  as  1 .2VCj  and 
do  not  ever  seem  to  slow  down  at  the  core  from  the  diffraction  event.  Note  the 
vortex  at  the  as  seen  by  the  by  the  streamlines  due  to  rapid  lateral  expansion. 
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Figure  37.  From  left  to  right  H2,  OH  mass  fractions,  pressure  (in  Pa) 

temperature  (in  K) 

The  following  two  figures  chronicle  the  detailed  re-ignition  of  the  front  for 
the  standard  case  with  time  sequenced  pressure  and  temperature  contours.  As 
the  wave  works  its  way  past  the  corner,  a  rarefaction  wave  originates  almost 
immediately  from  the  open  corner  begins  to  works  its  way  down  at  an  angle  of  30 
degrees  and  works  to  begin  a  decoupling  of  reaction  front  to  precursor  shock.  By 
frame  6,  the  rarefaction  wave  was  fully  into  the  reaction  zone  and  sweeps  across 
the  front  starting  to  flatten  it  and  adding  curvature. 

All  simulations  exhibited  a  decoupling  shock  and  reaction  front  from  the 
diffraction  front  followed  by  a  reinitiation  event  as  seen  in  the  above  Figures.  It 
should  be  noted  that  even  though  reinitiation  occurs,  the  shock  never  completely 
decoupled.  We  really  only  observed  localized  explosions  resulting  from 
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rarefaction  and  transverse  wave  collisions.  This  event  then  triggers  a  reaction 
across  the  unreacted  zone  or  shell.  Recoupling  of  the  shock  and  reaction  front 
then  spreads  over  the  entire  front.  Before  any  localized  explosion,  unburned 
reactants  slightly  affected  by  a  precursor  shock  find  themselves  in  a  spherical 
zone  between  a  flame  front  from  behind  and  a  transverse  wave  above.  In  this 
reinitiation  event,  it  appears  the  detonation  proceeds  transversely  across  the 
chocked  reactants  in  an  azimuthally  and  polar  manner  of  direction  completing  the 
reaction  and  reinitiating  the  entire  front,  as  seen  in  both  the  time  sequenced 
contours,  Figures  38  and  39. 
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Figure  38.  Time  sequence  pressure  contour  plot  with  streamlines  (in  Pa) 
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Figure  39.  Time  sequenced  temperature  contour  plots  (in  K) 
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VI.  EXPERIMENTAL  RESULTS 


A.  TESTING  OVERVIEW 

Table  5  depicts  an  overview  of  all  experiments  conducted  with  the  split 
flow  multi-stepped  configured  1 :3.5  aspect  ratio  channel.  The  tests  comprised  of 
detonating  hydrogen  and  air  at  atmospheric  conditions  for  various  diffraction  area 
ratios. 


Table  5.  Experimental  test  matrix 


Test# 

Test  Condition 

Equivalence 
Ratio  range 

Diagnostic 

1 

No  Step 

0.75-1.5 

Shadowgraph,  PIV 

2 

1  Step  1 :2.6  area  ratio 

0.75-1.5 

Shadowgraph 

3 

1  Step  1 :2.3  area  ratio 

0.75-1.5 

Shadowgraph 

4 

1  Step  1 :2.00  height  ratio 

0.75-1.5 

Shadowgraph, 

5 

2  Steps  1:2.3  &  1:2.3  area 
ratio 

0.75-1.5 

Shadowgraph, 

6 

2  Steps  1 :2  &  1:2  area  ratio 

0.75-1.5 

Shadowgraph,  natural 

emmision 

B.  PRESSURE  DATA 

Low  speed  pressure  measurements  were  primarily  used  to  calculate 
reactant  fuel  to  air  ratio  for  all  tests.  With  pressure  data,  equivalence  ratio  was 
calculated  based  on  the  choked  mass  flow  Equation  19,  with  main  air  flow  using 
a  0.235  inch  diameter  orifice  plate  while  a  0.157  inch  orifice  was  used  for  fuel 


=  Cr 


Ap, 

jV1 


f 

{  \ 

Vr 

,  P2 

1- 

Pi 

1 

VA  ) 

,  P,  j 

'Sc  ) 

2-y 

1 - 

{  R  J 

J- 

(19) 
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By  obtaining  a  ratio  of  mass  flow  rates  for  both  air  and  hydrogen,  the  equivalence 
ratio  could  be  determined  with  reasonable  errors  of  about  ±5%.  Discharge 
coefficient  Cd  for  both  orifice  plates  was  found  to  be  0.61  using  the  CRC 
handbook  for  mechanical  engineers.  These  results  were  also  in  good  agreement 
with  measurements  derived  PIV  flow  rates. 

High-speed  pressure  probe  data  was  acquired  using  Labview  and 
detonation  speeds  were  determined  by  measuring  the  inter  peak  distance  with 
time  taken  into  account  the  probe  spacing  of  0.241  m.  Typical  plots  obtained  can 
be  seen  in  Figure  40.  Typical  wave  speeds  ranged  from  1,800±50  m/s  to 
2,100±50  m/s  and  agreed  very  well  with  wave  speed  calculated  from  high  speed 
camera  imagery  and  CJ  wave  speeds  calculated  from  CEA. 


Figure  40.  High-speed  pressure  data 
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C.  RESULTS  OVERVIEW 

Figures  41  and  42  represents  a  summary  of  all  tests  performed.  The 
channel  performed  best  at  an  area  ratio  of  2  for  a  wide  range  of  equivalence 
ratios  for  both  single  and  double  step  configurations.  It  was  marginal  for  the  ratio 
of  2.3  and  did  not  transmit  for  any  area  ratios  of  2.6.  We  present  the  data  in 
terms  of  area  ratios  since  we  are  dealing  with  a  channel  that  is  not  perfectly  flush 
with  the  optical  section,  its  support  members  are  a  half-inch  thick  all  around  so 
there  is  some  diffraction  occurring  in  the  z-direction. 
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Figure  41 .  Single  Step  diffraction  results  with  Channel  Area  Ratio 

Detonations  transmit  into  the  double  step  area  ratio  of  2  for  three  of  the  six 
equivalence  ratios  tested  and  did  not  transmit  using  a  double  2.3  ratio,  seen  in 
Figure  42. 
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Double  Step  Diffraction  Results 
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Figure  42.  Double  step  diffraction  results  with  Channel  Area  Ratio 


Based  on  shadowgraph  imagery,  centerline  detonation  velocity  relative  to 

CJ  velocities  was  calculated  with  respect  to  distance  from  the  first  diffraction 

plane  for  a  number  of  equivalence  ratios  and  presented  on  6  of  the  following 

Figures.  Figures  43,  44  and  45  depict  wave  velocities  just  after  the  first 

diffraction  plane  at  three  equivalence  ratios  ranging  from  lean  to  rich  with  our  one 

step  configuration  for  varying  height  ratios.  The  area  ratio  of  2  performed  well 

under  all  conditions  as  observed  to  reinitiate  well  before  the  second  diffraction 

plane  (located  3.75  x/h  main  channel  heights  away)  coming  out  of  the  first 

diffraction  plane  as  wave  speed  were  close  to  CJ  velocity.  After  the  second  plane 

the  wave  diffracted  into  an  open  space  and  maintains  good  strength  for  another  2 

height  distances.  Area  ratios  of  2.3  did  not  fare  as  well,  as  the  reinitiation  seemed 

to  occur  at  a  later  distance  than  3.75  x/h.  Looking  at  the  0=1  case  there  was  a 

reinitiation  that  happened,  but  at  that  point  the  wave  did  not  have  enough  time  to 

completely  form  and  gain  enough  strength  to  survive  once  it  crossed  the  second 
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diffraction  plane.  It  seemed  that  a  good  step  height  to  length  ratio  was  2  times 
x/h,  so  long  as  the  wave  survives  the  initial  expansion  event.  For  our  area  ratio 
case  of  2.3,  the  second  step  height  is  2  inches  and  so  our  transition  length 
should  have  been  at  least  4  x/h,  allowing  enough  distance  for  a  detonation  to 
properly  reform.  Our  2.6  case  did  not  survive  the  first  diffraction  as  the 
expansion  waves  were  given  enough  time  and  distance  to  penetrate  the  entire 
front  enough  precluding  any  recovery,  in  all  reinitiation  cases. 
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Figure  43.  Detonation  velocity  vs.  run  distance  for  different  area  ratios 
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Figure  44.  Detonation  velocity  vs.  run  distance  for  different  area  ratios 
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Figures  46,  47  and  48  represent  the  results  for  both  double  step  setups  of 
Area  ratio  of  2  and  2.3.  The  double  area  ratio  of  2  successfully  transitioned  both 
steps  with  velocities  slightly  higher  in  the  3.75-6  x/h  zone  then  the  previous  case. 
The  addition  of  another  wall  to  created  another  shock  to  high  impedance  material 
interaction  creating  a  Mach  stem  reinitiation  zone  due  to  the  shock  reflection.  It 
should  be  noted  that  the  range  of  equivalence  ratios  for  survival  was  reduced. 
The  double  area  ratio  of  2.3  was  not  as  successful.  Longer  step  lengths  to 
accommodate  this  ratio  it  may  have  resulted  in  successful  reignition. 


1.2 

1.1 


1 


Double  Step  Diffraction  4>=0.8 

2nd  Diffaction 
plane 


1 


> 


o 

O) 

> 


0.9 

0.8 

0.7 


0.6 


0.5 


i  +  i 

+ 

*  V  * 

T 


0.4  — 

3.5  4  4.5  5  5.5  6 


+ 


6.5 


Channel  Distance  over  main  channel  height  x/h 


2 

2.3 


Figure  46.  Detonation  velocity  vs.  run  distance  for  different  area  ratios 
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Figure  47.  Detonation  velocity  vs.  run  distance  for  different  area  ratios 
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Figure  48.  Detonation  velocity  vs.  run  distance  for  different  area  ratios 
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D.  DIFFRACTION  PREDICTION  MODEL 

A  simple  diffraction  failure  model  was  tested  and  validated  with  these 
experimental  results.  Based  on  previous  results  both  detonation  cell  size  for 
hydrogen  and  air  with  respect  to  equivalence  ratio  and  detonation  velocities 
could  be  predicted  for  the  various  reactants  as  seen  in  Table  6.  Knowing  these 
parameters,  our  detonation  channel's  geometric  dimensions  and  how  critical  our 
channel  was  for  an  aspect  ratio  of  3.5  we  could  predict  the  velocity  deficit  under 
these  conditions  [24],  We  could  also  calculate  either  rarefaction  angle  from 
disturbance  to  undisturbed  detonation  a  or  a  reaction  zone  velocity  Vrz  to  the  first 
step  from  Skews  method  [25],  depending  on  what  data  was  available.  In  our 
case  we  calculated  Vrz  using  Equation  20  and  measured  a  based  on  the  Shultz's 
high  speed  shadowgraph  images  of  Hydrogen  and  air  and  our  own  shadowgraph 
images  which  was  found  to  be  15  degrees  [8].  Knowing  our  step  distance  and 
the  distance  required  for  total  failure  of  the  front,  we  could  predict  the  health  of 
the  front  by  the  time  the  wave  hits  the  step  wall.  Degree  of  criticality  in  Table  6 
refers  to  how  many  detonation  cell  widths  we  were  across  the  cross  section  in  a 
particular  case  with  respect  to  the  critical  height  for  a  3.5  aspect  ratio  channel 
which  is  7A.  For  most  cases,  we  were  about  5  times  smaller  than  this  so  called 
critical  length.  With  this  information,  we  could  apply  a  velocity  deficit  linked  to 
criticality  [7]. 
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Figure  49.  Schematic  of  a  diffracted  shock  (From  [8]) 


Based  on  our  experiments,  it  appeared  that  once  the  entire  front  had 
curvature  the  diffracted  chock  had  no  chance  of  survival.  However,  if  some 
portion  of  the  reaction  zone  was  still  coupled  to  the  pre-cursor  shock  there  was  a 
chance  of  survival  if  a  Mach  stem  developed  quickly  enough  from  the  wall  area 
and  had  time  to  progress  into  the  rest  of  the  front. 

Our  case  with  an  area  ratio  of  2.6  failed  every  time  experimentally  and 
was  predicted  in  our  model.  We  also  find  experimentally  that  the  front  has 
already  been  completely  disturbed,  seen  in  Table  7.  Whereas  area  ratio  of  2 
detonations  survived  every  single  step  expansion  experimentally  and  we  predict 
that  this  wave  was  about  80%  failed. 
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Table  6. 


Predicted  numerical  results  with  a=15°,  Acr=7 


<p 

A  (in 
cm) 

A's  across 

area 

Acr/A's 

across 

(criticality) 

area 

ratio 

(AR) 

criticality 

/AR 

0.8 

2 

1.25 

5.60 

2.00 

2.80 

0.8 

2 

1.25 

5.60 

2.30 

2.43 

0.8 

2 

1.25 

5.60 

2.60 

2.15 

1 

1.5 

1.67 

4.20 

2.00 

2.10 

1 

1.5 

1.67 

4.20 

2.30 

1.83 

1 

1.5 

1.67 

4.20 

2.60 

1.62 

1.15 

1.7 

1.47 

4.76 

2.00 

2.38 

1.15 

1.7 

1.47 

4.76 

2.30 

2.07 

1.15 

1.7 

1.47 

4.76 

2.60 

1.83 

Table  7.  Predicted  numerical  results  cont. 


<P 

area 

Vcj 

V  deficit 

V 

Time  t  to 

% 

Exp 

ratio 

(m/s) 

(m/s) 

(m/s) 

wall  (s) 

failed 

results 

0.8 

2.00 

1866 

1810.02 

905.01 

4.88E-05 

84.19 

passed 

0.8 

2.30 

1866 

1810.02 

905.01 

5.61  E-05 

96.76 

passed 

0.8 

2.60 

1866 

1810.02 

905.01 

6.35E-05 

109.52 

failed 

1 

2.00 

1966 

1946.34 

973.17 

4. 54  E-05 

84.19 

passed 

1 

2.30 

1966 

1946.34 

973.17 

5.22E-05 

96.76 

passed 

1 

2.60 

1966 

1946.34 

973.17 

5.91  E-05 

109.52 

failed 

1.15 

2.00 

2015 

1974.7 

987.35 

4.48E-05 

84.19 

passed 

1.15 

2.30 

2015 

1974.7 

987.35 

5.15E-05 

96.76 

passed 

1.15 

2.60 

2015 

1974.7 

987.35 

5.82E-05 

109.52 

failed 

E.  TEST  CONDITION  1 :  COLD  FLOW  AND  WIDE  OPEN  DIFFRACTION 

Test  condition  1  included  all  experiments  that  would  help  understand  all 
the  conditions  the  detonation  front  would  be  exposed  to  in  the  simplest 
configuration,  i.e.,  diffraction  into  a  wide  opening  without  the  stepped  sections  in 
the  optical  section.  Tests  included  basic  PIV  tests  of  the  cold  flow  and 
diffractions  into  an  open  channel  as  well  shadowgraph  shots. 
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Cold  flow  tests  were  conducted  to  better  understand  the  pre  detonation 
flow  conditions  and  seeding  quality  of  the  flow.  Figure  49  was  a  typical  test  run 
with  main  air  flow  mass  flow  rate  of  0.017  kg/s  depicting  camera  bit  count 
(intensity  level)  of  the  seed.  The  laser  sheet  illuminates  an  80  mm  by  80  mm 
zone  in  the  flow  immediately  following  the  diffraction  plane,  flow  moves  from  left 
to  right.  Note  that  these  counts  were  slightly  high  as  the  bulk  of  the  seed  was 
biased  towards  lower  main  detonation  channel.  Although  seed  density  was 
higher  in  this  region  than  required  at  these  low  speed  fills,  it  was  found  that  such 
levels  were  required  for  sufficient  seed  density  required  for  good  vector 
correlation  for  wave  speeds  in  the  order  of  2,000  m/s  as  such  waves  would  tend 
to  disperse  the  seed  after  the  front,  as  particles  can  be  accelerated  to  as  much 
as  1 ,000  m/s  just  after  the  front  applying  significant  seed  density  variations. 
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Figure  50.  Seed  bit  counts  with  low  speed  cold  flow  run 


The  postprocessed  vector  plot  for  Figure  50  can  be  seen  at  Figure  51 . 
Flow  on  average  was  very  turbulent  and  fluctuated  significantly  during  the  filling 
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process.  On  average,  velocities  in  the  top  channel  were  higher  by  a  factor  of  2 
compared  to  flow  in  the  lower  channel.  A  shear  layer  of  lower  velocity  can  be 
found  in  between  both  channels. 
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Figure  51 .  Cold  Flow  PIV  velocity  vector  plot 


The  following  Figures  depict  PIV  vector  plots  of  one  snapshot  of  the 
detonation  front  without  any  steps  with  a  0=1 .  The  first  Figure  is  frame  two  of 
the  two  frame  camera  showing  the  good  seed  density  in  and  around  the 
diffracting  detonation  front.  A  pressure  of  5  psig  was  applied  to  the  seeding 
system  for  reasonable  seeding  density.  Inter  laser  time  was  calculated  in  Davis, 
Lavision's  in  house  PIV  code  based  on  estimated  800  m/s  particle  velocity  behind 
the  front,  resulting  in  a  710  ns  time  deference  between  frames  dt.  The  trigger 
delay  was  set  to  230  ps  and  triggering  was  off  of  the  first  Kistler  pressure  probe 
with  a  laser  power  of  50%  of  max. 
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velocity  m/s 


Figure  52  was  the  post  processed  vector  plot  of  the  flow  field.  We 
observed  particle  speeds  in  the  800  m/s  range  near  the  core  dropping  to  as  low 
as  400  m/s  if  we  follow  around  the  perimeter.  Unfortunately,  with  such  a  small  dt 
we  missed  capturing  lower  end  velocities  and  similar  test  runs  would  be  needed 
with  various  dt's  to  get  the  whole  event.  Additionally,  because  we  used  a  low 
frame  rate  PIV  speed  system  to  capture  detonations  at  various  diffraction  times 
many  runs  would  have  been  required. 
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Figure  52.  Seed  bit  counts  with  high  speed  settings 


The  postprocessed  plot  was  quite  noisy  around  the  front.  It  is  suspected 
that  fine  tuning  camera  focus  and  seed  density  would  have  improved  the 
proportion  of  noise  and  we  would  have  obtained  better  correlation  with  calculated 
vectors. 
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Figure  53.  Post  processed  vector  PIV  vector  plot  of  a  diffracting  detonation 

front 


Shadowgraph  sequences  were  also  taken  without  any  step  sections  this 
can  be  seen  in  Figure  54.  All  images  were  shot  at  150,000  frames  per  second 
with  a  total  of  17  frames  with  the  Ultra  17  high-speed  camera  system. 
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Figure  54.  Shadowgraph  sequence  of  diffraction  into  unconfined  area 
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F. 


TEST  CONDITION  2:  SINGLE  STEP  AREA  RATIO  2  WITH 
SHADOWGRAPH 


Figure  55.  One  step,  0=1 , 2.00  A2/Ai 


Figure  56.  Frame  4  From  above,  x/h  1-3.5 
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G. 


TEST  CONDITION  3:  SINGLE  STEP  AREA  RATIO  2.3  WITH 
SHADOWGRAPH 


Figure  57.  One  step,  0=1 , 2.30  A2/Ai 
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Figure  58.  Frame  8  from  above,  x/h  2-4.5 
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H. 


TEST  CONDITION  4:  SINGLE  STEP  AREA  RATIO  2.6  WITH 
SHADOWGRAPH 


Figure  59.  One  step,  <t>=1 , 2.6  A2/A-i 
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I. 


TEST  CONDITION  5:  DOUBLE  STEP  AREA  RATIO  2.00  WITH 
NATURAL  EMISSION 


The  following  two  images  were  taken  with  high-speed  camera  looking 
directly  at  the  optical  section  without  laser  illumination.  A  reinitiation  event  can 
be  observed  moving  azimuthally  across  what  was  once  a  decoupled  reaction 
zone.  In  Figure  62,  the  front  is  still  planar  at  a  distance  of  5.5  x/h. 


Figure  61.  Natural  emission,  0=1, 2.00  A2/A1,  x/h  1  -  6 


Figure  62.  Natural  emission,  0=1, 2.00  A2/Ai,  x/h  1  -  6 
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J. 


TEST  CONDITION  6:  DOUBLE  STEP  AREA  RATIO  2.3  WITH 
SHADOWGRAPH 


Figure  63.  Two  step,  0=1 , 2  A2/A-i 


Figure  64.  Frame  12  from  above 
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VII.  FUTURE  WORK 


Future  work  could  encompass  both  computational  and  experimental 
aspects.  Computationally,  the  next  step  would  be  to  rerun  these  results  with  a 
more  refined  structured  mesh  and  to  develop  a  better  chemistry  set  for 
detonations. 

Experimentally,  longer  step  sections  could  be  manufactured  to  ensure 
proper  full  transition  during  diffractions  events.  Also,  different  fuels  could  be 
used  to  test  the  predictions  model  as  well  as  a  performance  comparison  to 
hydrogen.  Additional  PIV  work  is  required  once  Lavision  rebuilds  the  lab's 
existing  camera.  With  the  new  optical  section  more  diffraction  experiments  could 
be  conducted  using  CH*  and  PLIF. 
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VIII.  CONCLUSIONS 


A  new  optically-accessible  test  section  was  designed  and  built  to  image 
the  detonation  diffraction  of  hydrogen  and  air  mixtures  across  multi-step 
configurations.  Shadowgraph  and  PIV  imaging  techniques  were  used  to 
determine  shockwave  profiles,  velocities,  and  re-ignition  mechanisms  throughout 
the  diffraction  process.  Computational  models  were  investigated  but  failed  to 
produce  any  useful  information  due  to  the  artificial  high  reactivity  of  the  mixtures 
for  the  various  reduced  chemistry  kinetics  and  also  inadequate  mesh  resolution 
utilized. 

Experiments  revealed  that  hydrogen/air  detonations  could  survive  single 
and  double  diffraction  across  an  area  ratio  of  2.0  for  an  equivalence  ratio  range 
of  0.8  to  1.14.  Successful  detonation  diffraction  was  observed  for  area  ratios  of 
2.3  but  only  for  the  primary  diffraction  event  and  only  for  equivalence  ratios  near 
stoichiometric  for  the  given  step  length.  Increasing  the  step  length  would 
improve  the  likelihood  of  subsequent  diffraction  events  by  allowing  the  detonation 
to  reform  before  encountering  another  diffraction  condition.  No  successful 
detonation  diffraction  reinitiations  were  observed  for  area  ratios  of  2.6. 
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APPENDIX  A.  NASA  CHEMICAL  EQUILIBRIUM  OUTPUT 


NASA-G LENN  CHEMICAL  EQUILIBRIUM  PROGRAM  CEA2  ,  FEBRUARY  5,  2004 

BY  BONNIE  MCBRIDE  AND  SANFORD  GORDON 
REFS:  NASA  RP-1311,  PART  I,  1554  AND  NASA  RP-1311,  PART  II,  1936 


A  *  -A  A  A  A  A  -J.  *  A-  A  -A  A  -A  A  A  A  A  A  A  A  A  Ai  A  A  A-  i-  A  A  i  A  A  -A  AA.A-J.AA  *******  **********  A  -A  A  A  -A  A  A  A  A  A  A  A  A  A 
*******  * 


output  massf 

prob  case=00Q02142  det 

phi— 1 

p (atm) -1 
output  short 
output  tra ce=le-I 0 

reac 

oxid  Air  t,k=  258  wt*=  100. 
fuel  H2  t , k=  258  wt%=  100. 
end 

DETONATION  PROPERTIES  OF  AN  IDEAL  REACTING  GAS 

CASE  =  00002142 


TEMP 


REACTANT 


K 

OXIDANT  Air 

298.000 

FUEL  H2 

298 . 000 


VFT  FRACTION 
(SEE  NOTE) 

1 . 0000000 
1.0000000 


ENERGY 


KJ/KG-MCL 


-125 .855 


-4.325 


Q/F=  34.29623  %FUEL=  2.833164  RrEQ.RATIO=  1.000000  PHI ,  EQ  .  RATIO— 

1.000000 


UNBURNED  GAS 


PI,  BAR  1.0132 
Tl,  K  238.00 
HI,  KJ/KG  -4.36 
Ml,  ( 1/n)  21.008 
GAMMA  1  1.4015 
SON  VEL1, M/SEC  406.6 


BURNED  GAS 
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T,  K 

2543.63 

RHO,  KG/ CU  M 

1.5500  0 

H,  KJ/KG 

1333.33 

U,  KJ/KG 

313.50 

G,  KJ/KG 

-25333, 7 

S,  KJ/  (KG)  (K) 

10 . 5G80 

M,  ( 1/n) 

24.008 

(dLV/dLP) t 

-1.00554 

( dL V / dLT ) p 

1.2062 

Cp,  KJ/  (KG)  £K) 

3.3542 

GAMMA  s 

1 . 1637 

SON  VEL,  M/SEC 

1005.2 

DETONATION  PARAMETERS 

P/Pl 

15 . 555 

T/Tl 

3.078 

M/ML 

1 . 1428 

RHO/ RHO 1 

1 . 0042 

DET  MACH  NUMBER 

4 . 8335 

DET  VEL, M/SEC 

1565.1 

MASS  FRACTIONS 

*  A_r 

1 . 2550-2 

*CO 

1 . 3053-4 

*C02 

2. 6525-4 

CO  OH 

2 . 6145-5 

*H 

2 . 4760-4 

HCO 

1 . 0225-9 

HNCO 

2 .3  95-10 

HNO 

3 . 1235-6 

HN02 

8 . 9835-7 

HND3 

1 . 020-10 

HO  2 

1 . 5250-5 

AH2 

2 . 6254-3 

HCOOH 

3 . 735-10 

H20 

2 . 2035-1 

H2D2 

2 . 3000-6 

*N 

1 . 1146-6 

+  NH 

4 . 5040-7 

NH2 

2 . 7000-7 

NH3 

3 .4965-7 

NH20H 

4 . 217-10 

*  NO 

5 . 4235-3 

NO  2 

6 . 6335-6 

*N2 

7. 2938-1 

N2  0 

2 . 7711-6 

*0 

1 . 3571-3 

*OH 

1 . 3522-2 

+  02 

1 . 0117-2 

03 

6 . 2070—5 
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APPENDIX  B.  CAD  DRAWINGS 


The  following  set  of  Figures  were  all  generated  with  the  Solidworks  2010 
mechanical  modeling  software  package.  Optical  section  and  the  new  head  end 
inlet  flange  were  made  from  Type  303  Stainless  Steel  in  order  to  meet  the  rigidity 
and  strength  requirements  for  the  engine  test  section  while  also  being  the  most 
readily  mechineable?  for  an  austenitic  stainless  steel.  All  internal  components 
were  made  from  6061  Al.  All  units  are  in  inches  and  degrees.  Tolerances  for 
linear  dimensions  are  ±0.01"  and  radial  dimensions  are  ±0.05°. 
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Detonation  Channel 


Figure  65.  Overall  detonation  channel  assembly 


Figure  66.  Exploded  view  of  Optical  assembly 
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Figure  67.  Back  plate  (1  of  2) 
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8.500 


Figure  68.  Back  plate  (2  of  2) 
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36X  0  .28  THRU  ALL 


32  X  0  .21  \F  .65 
1  /4-28  UNF  T?  .50 


Figure  69.  Outer  wall  optical  left  and  Right  (1  of  4) 


22X  0  .21  ^  .65 
1/4-28  UNF  ^  .50 


90 


DETAIL  D 
SCALE  1  :  1 


nyuit;  iv.  wuitM  wciii  upuucii  i_eu  anu  myiu  ui  h) 
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Figure  7 1 .  Outer  Wall  optical  Left  and  Right  (3  of  4) 
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Figure  72.  Outer  Wall  optical  Left  and  Right  (4  of  4) 
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Figure  74.  Inner  wall  optical  (2  of  2) 


95 


96 


__  i  I  on  __ 

m  1  1  ->-AJ  ■■ 

.45 

—  n  C\i=: 

f  .25  - 

.4A 

-i 

■  ^  o  y  n  — . 

r  y-/u  *i. 

Figure  76.  Glass  Holder  optical  (2  of  2) 
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Note  1 


Note  1 :  optical  polished  surface  or 
both  front  and  back 

Note  2:  all  sharp  edges  and  corner: 
removed 


Figure  78.  BK7  glass  for  main  optical  section 


Figure  79.  Bottom  plate  assembly  for  PIV  laser  sheet 
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Figure  80.  Bottom  plate  with  BK7  glass  slot  PIV  laser  sheet  (1  of  2) 
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SCALE  1  : 2 
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4x  1/4-20  0  0.75  nF 


E1 IMENSIO  NS  ARE 

TOLERANCES: 

FRACTIONAL:; 


Rocket  Lab 


Figure  84.  Aluminum  Insert 
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Insert  Front 


Front  Adapter 


Figure  85.  Split  flow  insert  assembly 
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Cl 


Figure  88.  Insert  top  (1  of  2) 
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DETAIL  C 
SCALE  1  : 4 


DETAIL  A 
SCALE  1  :2 


1  \ 

0.26 

.50 

"-“-B 

i|.t 

*tl - 

- -fc  -2- -=_1_ -=_■=_■= 

1“ 

3*45 

-i  1  ^  tr 

^  1  Jz.  ^ 

0.40 

Figure  90.  Insert  bottom  (1  of  2) 
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DETAIL  6 

SCALE  1  *2 

Figure  91 .  Insert  bottom  (2  of  2) 
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Figure  92.  Insert  left  (1  of  2)  (mirrored  for  right  side) 
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Figure  93.  Insert  left  (2  of  2)  (mirrored  for  right  side) 


Figure  94.  Insert  front 
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SECTION  D-D 
SCALE  1  ;  2 


Figure  96.  Insert  ramp  holder  (1  of  2) 
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Figure  97.  Insert  ramp  holder  (1  of  2) 
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JO-54) 


00.1 1  .  0.10 


Small  ramps 
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Figure  99.  Step  insert  sample  (steps  made  for  1:1.25,  1:1.5  and  1:1.75) 
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6x0  0.280  ^  0.253 
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1 

.50 
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I 

1.32 

r 

t 

- -  04.50  - - 

ingiter  plug  tap 
1 1  /I 6-24 
0.25  deep 


inner  dia  0.56 
thru 


Figure  1 02.  Inlet  manifold  (1  of  2) 
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APPENDIX  C.  STANDARD  OPERATING  PROCEDURE 


Test  Cell  #1  Standard  Operating  Procedures  (SOP) 

Combustor  Start  Up 
(Modification  Date  4  Nov  10) 

INITIAL  SAFETY  CHECKS 


1 .  Notify  all  lab  personnel  of  intention  to  make  test  cell  1  live. 

2.  Turn  ON  control  console 

3.  Turn  ON  warning  lights 

4.  Cell  #1  EMERGENCY  SHUTDOWN  BUTTON  (Control  Room)- VERIFY 
IN 

INITIAL  PREPARATIONS 

1.  Notify  the  Golf  Course  (x21 67)  (Only  required  if  Hot  Fire  Test  is 
conducted) 

2.  Open  Test  Cell  Door 

3.  Igniter  Control  (Test  Cell)-VERIFY  OFF  (Red  Button  Out) 

4.  PXI-1000B  Rack  (Test  Cell)-VERIFY  ON 

5.  Shop  Air-VERIFY  >  100  PSI 

6.  Shop  Air  Valve  (Test  Cell)-VERIFY  OPEN 

7.  1 15  VAC  Control/Cell  #1  Switch  (Control  Room)-ON 

8.  28VDC  Power  Supply/Cell  #1  Switch  (Control  Room)-ON 

9.  Open  LABVIEW  file  "TCI  ControllerMike" 

TESTING  SET-UP 

1 .  Inspect  optical  section  and  mirrors  for  cleanliness 

2.  Commence  Argon-Ion  laser  start-up  procedure  (see  Arg-lon  SOP) 

3.  Commence  PIV  NG-Yag  double  pulse  laser  start  up  (PIV  SOP) 

4.  PIV  software  controller  -  ON 

5.  Kistler  Amplifier  Power-ENSURE  OFF 

6.  Kistler  Leads-CONNECT 

7.  Exhaust  Tube-VERIFY  PROPER  POSITION 

8.  Notify  all  personnel  that  gasses  and  TESCOM  will  be  enabled. 

9.  Test  Cell  #2  and  #3  Node  1  Air  Isolation  Valve  (Test  Cell  #2)-VERIFY 
CLOSED 

**NOTE:  This  valve  maybe  left  open  only  if  Test  Cell  #2  or  #3  is  configured 

for  active  testing** 
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10. TRANSDUCER  and  TESCOM  Power  Switch  (Test  Cell  #2)-ON 
1 1  .Set  0  (Zero)  pressures  on  ER3000  (Control  Room)  for  the  following: 

a.  Node  3  (Hydrogen) 

b.  Node  1  (Main  Air) 

c.  Node4  (Seeder  Air) 

12.  Main  HP  Air  Jamesbury  Valve  (Outside  Test  Cell)-OPEN  SLOWLY 

13.  Power  Strip  (above  PXI-1000B)-VERIFY  ON 

14.  Igniter  Control  Light  5V  logic  (Red  LED  CRYDOM  2nd  Column,  2nd  down 
in  PX-1000B  Rack)-VERIFY  OUT 

**DANGER:  IF  RED  LIGHT  IS  ENERGIZED,  MUST  CLEAR  USING  LAB  VIEW 
BEFORE  CONTINUING  TO  PREVENT  PREMATURE  IGNITION** 

15.  Igniter  Control  (Test  Cell)-PUSH  RED  BUTTON  IN 

16.  Igniter  Control  Startup  Diagnostics-OBSERVE  COMPLETION  OF 
DIAGNOSTICS  (Verify  energy  level  setting  reads  2.03  J) 

17.  Main  HP  Air  Isolation  Valves  (2)  (Located  in  Test  Cell)-OPEN  SLOWLY 

**DANGER:  OPEN  VALVES  SLOWLY  TO  PREVENT  RAPID 
PRESSURIZATION  OF  DOWNSTREAM  LINES** 

18.  Node  3  (Hydrogen)  Shop  Air  Valve  (Above  Bottle  in  Bottle  Room)-OPEN 

19.  Hydrogen  Bottle  Isolation  Valve  (On  Bottle)-OPEN  SLOWLY 

**  VERIFY  ADEQUATE  ETHYLENE  PRESSURE  FOR  TESTING  ON 
DOWNSTREAM  GAGE  LOCATED  IN  BOTTLE  ROOM** 

PRE-TESTING 

1 .  SA-5  Camera-TURN  ON  AND  REMOVE  COVER 

2.  Start  SA-5  Software 

3.  Verify  Image  on  SA-5  Software 

4.  PIV  Camera  -  TURN  ON  AND  REMOVE  COVER 

5.  Start  PIV  software 

6.  Verify  Desired  Trigger  Types  and  Valve  Durations 

7.  Determine  Desired  Fuel  and  Air  Pressures  (Mass  Flow  Choke 

Calibrations.xls) 

8.  Set  Required  Pressures  on  ER3000 

a.  Node  3  (Hydrogen) 

b.  Node  1  (Main  Air) 

c.  Node  4  (Seeder  Air) 

9.  Ready  Digital  Camera 

10.  Switch  Monitor  to  Test  Cell  #1 
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TESTING 


1 .  Clear  All  Test  Cells  and  Verify  with  Head  Count 

2.  Flashing  Red  Light  and  Siren-ENERGIZE 

3.  Verify  Golf  Course  is  CLEAR 

4.  Digital  Camera  RECORD 

5.  Camera-ARM 

6.  In  LABVIEW  Enable  Facility  Button-ON 

7.  Test  Cell  #1  Emergency  Shutdown  Button-TURN  CLOCKWISE 

8.  In  LABVIEW  Start  Button-CLICK  TO  START 

WHEN  TESTING  COMPLETE 

9.  Set  Node  1  Pressure  to  0  (Zero) 

10.  In  LABVIEW  Turn  Off  Button-CLICK  TO  SECURE 

1 1 .  In  LABVIEW  Enable  Facility  Button-VERIFY  OFF 

12.  Test  Cell  #1  Emergency  Shutdown  Button-PUSH  IN 
13.Siren-OFF 

14.  Digital  Camera-Stop/Pause 

15.  Save  SA-5  Image 

16.  Save  PIV  data 

**NOTE:  If  Further  Testing  is  Required,  re-perform  steps  1-15  of  the 

Testing  Section** 

POST-TESTING 

1 .  Close  fuel  isolation  valves,  rerun  labview  code  to  cycle  fuel  valve 

2.  Igniter  Control  (Test  Cell)-PUSH  RED  BUTTON  OUT 

3.  Flow  additional  Main  Air  to  purge  any  remaining  fuel 

4.  Set  ER300  Node  3  (Ethylene)  to  0  (Zero) 

5.  Verify  pressures  on  ER3000  are  set  to  0  (Zero)  on  the  following: 

a.  Node  3  (Hydrogen) 

b.  Node  1  (Main  Air) 

c.  Node  4  (Seeder  Air) 

6.  Main  HP  Air  Isolation  Valves  (2)  (Located  in  Test  Cell)-CLOSE 

7.  Main  HP  Air  Jamesbury  Valve  (Outside  Test  Cell)-CLOSE 

**NOTE:  If  Further  Testing  will  be  accomplished  with  a  different  Ramp 
obstacle  configuration,  return  to  the  TESTING  SET-UP  Section.  If  not 

continue  to  step  7.** 

8.  SA-5  Camera-TURN  OFF  AND  INSTALL  COVER 

9.  PIV  Camera  -  TURN  OFF  AND  INSTALL  COVER 
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10.  TRANSDUCER  and  TESCOM  Power  Switch  (Test  Cell  #2)-OFF 

1 1 .  Close  Test  Cell  Door 

12.  Node  3  (Hydrogen)  Shop  Air  Valve  (Above  Bottle  in  Bottle  Room)-CLOSE 

13.  Secure  Bottle  Room 

14.  Exit  out  of  SA-5  and  PIV  software 

15.  EXIT  out  of  LABVIEW 

10.28VDC  Power  Supply/Cell  #1  Switch  (Control  Room)-OFF 
1 1 . 1 15  VAC  Control/Cell  #1  Switch  (Control  Room)-OFF 
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